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DURCH KRISTALLE VERURSACHTE 
ZERSTORUNGSWIRKUNGEN IN POROSEN 
MATERIALIEN 
Von 
Z. GYULAI 


INSTITUT FUR EXPERIMENTELLE PHYSIK DER TECHNISCHEN UNIVERSITAT 
FUR DIE BAUINDUSTRIE, BUDAPEST 


(Eingegangen: 25. VI. 1954) 


An porésen Kérpern — so an gewissen Porzellanen, Fayencen, keramischen Materi- 
alien — wachsen nach Aufsaugung der NaCl-Lésung nadelférmige und wollartige NaCl- 


- Kristalle. Steht der Probekérper sehr lange Zeit (mehr als 2 Jahre) an freier Luft, so entstehen 


an gewissen Stellen Aufquellungen, Abbréckelungen und Spriinge. Diese Zerstérung des Probe- 
kérpers kann durch die Entstehung von Kristallen erklart werden, die in den Porenraumen 
zu wachsen beginnen. Es wird auch auf die Bedeutung dieser Erscheinung bei der zeitlichen 
Verwitterung der Baumaterialien (Ziegel, Steine, Zemente) hingewiesen. 


1. Am I. Kongress der Ungarischen Physiker im August 1953 berichtete 
der Verfasser ttber die kiinstliche Ziichtung von nadelférmigen und wollartigen 
NaCl-Kristallen. Diese Ziichtungsversuche sind seit mehreren Jahren im Gange, 
und es wurden verschiedene Materialien ausprobiert, ehe es gelang, die heutige 
sog. Massenproduktion zu erreichen. Der Ausdruck »Massenproduktion« wird 
hier deshalb gebraucht, weil es sehr schwierig ist, mit Nadelkristallen zu arbeiten, 
deren Dickenabmessungen in einen Bereich von 1—10 yw fallen, und weil geniigend 
Kristallmaterial zur Verfiigung stehen muss, wenn man vermeiden will, dass 
die wahrend des Versuches auftretenden Briiche usw. eine Verzégerung der 
Arbeit hervorrufen. 

Im nachstehenden soll — gleichsam als Nebenergebnis — itiber gewisse 
Zerstérungs- bzw. Zerbriéckelungserscheinungen berichtet werden, die auf dem 
als Grundmaterial dienenden porésen Stoff auftreten. Diese Beobachtungen 
bilden auch eine Erweiterung des auf dem Kongress gehaltenen Berichtes in 
eine andere Richtung. Es wurde nimlich angenommen, dass das Kristallwachs- 
tum nur an der Oberflache eines porésen Stoffes vor sich gehen kann, da die 
Lésung dort zu verdunsten vermag. EKingehende Beobachtungen ergaben jedoch, 
dass der Verdunstungsprozess etwas komplizierter ist, als urspriinglich ange- 
nommen wurde. 

2. Die bei diesen Versuchen gemachten Beobachtungen sind folgende : 
es wurde in einen Fayenceteller gesattigte Kochsalzlésung gegossen und unge- 
fahr 24 Stunden im Teller stehen gelassen, worauf der Teller mit einem Tuch 
getrocknet wurde. Auf der Innenseite des Tellers begannen nach einigen Monaten 
die normalen wollartigen Kristalle zu wachsen. Nach zwei Jahren — wahrend 
denen der Teller standig in trockener Zimmerluft gehalten wurde —. zeigten 
sich jedoch an seiner Oberfliche gréssere Erhabenheiten bzw. Aufquellungen. 


T Acta Physica IV/3. 
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Ein solcher Fall ist in Abb. 1 zu sehen. Abb. 1/a zeigt den ganzen Teller, am 
Rande an einigen Stellen mit den Kristallbiindeln, die selbst nach zeitweiligem 
Abbrechen, seit zwei Jahren ein stindiges Wachstum aufweisen. Am Boden des 
Tellers ist eine Flache von der Grésse einer 1-Forint-Minze aufgequollen, 


b) 


Abb. 1. Fayenceteller mit aufgequollenem Teil. a) Rand des Tellers mit Kristallhaufen, 
in der Mitte ein aufgequollener Teil. b) Aufgequollener Teil in starkerer Vergrésserung, 
neben der gréssten Erhebung noch ein kleinerer Teil ausgebrochen. 


daneben sind an zwei kleineren Teilen und auch an anderen Stellen des Tellers 
ausgesprochene Erhebungen zu sehen. Am Teller ist auch ein Sprung entstanden, 
der erst in der letzten Zeit festgestellt werden konnte. 

Diese Sprung- bzw. Spalterscheinungen dirften auf die sich im Inneren 
der Fayence bildenden Kristalle zuriickzufiihren sein. 

Auf welche Weise kénnen im Inneren der Fayence Kristalle wachsen? 
Hier sind allerkleinste Offnungen vorhanden, deren Abmessungen wahrschein- 
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lich in einer Gréssenordnung von 1 u bis 0,1 « und vielleicht sogar bis 0,01 u 
liegen. Wenn eine als kreisformig angenommene Offnung 0,01 gross ist, 
d.h. einen Durchmesser von 100 A aufweist, dann steht den Wassermolekiilen 
noch immer reichlich Raum zur Verfiigung. Setzt man also an der Wand der 
inneren Poren eine adsorbierte Fliissigkeitsschicht von der Dicke eines oder 
mehrerer Molekiile voraus, so ist iiber dieser Schicht ein freier Raum von etwa 
100 Molekiile Héhe zur Verdunstung der Wassermolekiile vorhanden. Dieser 
Raum fillt sich mit gesattigtem Wasserdampf, wobei die Salzionen an der 
Seitenwand der Poren eine Flachenwanderung ausfiihren kénnen. Obwohl 
besonders zu betonen ist — was auch bereits im Vortrag auf dem Kongress 
kurz gestreift wurde —, dass die Angaben der adsorbierten Schicht (Ionen- 
-konzentration, Dampfdruck) andere sind als die fiir die gesattigte NaCl-Liésung 
geltenden betreffenden Angaben, so unterliegt es keinem Zweifel, dass sich 
bei Ubersattigung der adsorbierten Schicht die Ionen aus ihr in Form von 
Kristallen ausscheiden. Es besteht somit keinerlei Hindernis, dass Kristalle 
auch in solchen Poren wachsen. 

Dieser Prozess nimmt natiirlich eine lange Zeit vom Versuchsbeginn an 
gerechnet, d.h. vom Bedecken der Fayence mit Salzlésung, in Anspruch. 
Wahrscheinlich fillt sich am Anfang ein Teil der Poren vollstandig mit Lésung 
an, wahrend in einem anderen Teil von ihnen Lufteinschliisse zuriickbleiben, 
weil gewisse Teile infolge der in die iibrigen Poren auf Grund der Kapillaritat 
eindringenden Fliissigkeit abgesperrt werden. 

Nach Ansetzen des Versuches wurde — wie bereits erwahnt — der Fayence- 
teller mit einem Tuch abgetrocknet. Von diesem Zeitpunkt an beginnt in den 
Poren eine iiberaus verwickelte Bewegung der Fliissigkeit und ihres Dampfes. 
Es ist offenbar, dass sich infolge der Oberflachenwanderung auch an den Wanden 
der leeren Poren eine adsorbierte Schicht ausbildet. Infolge der Temperatur- 
anderung treten Volumveradnderungen auf, was gleichfalls die Bewegung der 
Flissigkeit und des Dampfes in den Poren férdert. Die Temperaturanderungen 
bewirken eine Veranderung der Dampfspannung und zugleich damit auch der 
Dampfdichte: es setzt also ein Spiel von Verdunstung und Niederschlag ein. 
Auch die Veranderung des Aussenluftdruckes spielt eine Rolle, indem sie eben- 
falls die Bewegung der Fliissigkeit und des Dampfes auslést. Alle diese Prozesse 
haben zur Folge, dass sich der Ort und die Konzentration der Gesamtheit 
der Ionen verandern. 

Nach all dem ist es also offenkundig, dass die Wassermolekile an den 
freien Enden der Poren — gleichgiiltig ob sich der Fayence-Probekérper in 
einer feuchten oder trockenen Umgebung befindet* — aus dem Porensystem 
_austreten, wahrend die zu ihnen gehérenden Ionen dort bleiben. In den bisherigen 
Versuchen wurden nur jene Ionen untersucht, die sich an der Oberfliche des 


* Ausfihrlicher s.: Z. Gyulai: Magyar Fizikai Folydirat, 2. 1954, 371—391. 
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Probekérpers in Form von nadelférmigen, geraden oder wollartigen, krummen 
Kristallen ausgeschieden hatten. 

Es ist offenbar, dass diese an der Oberfliche ausgeschiedene Substanz 
die Mehrzahl der vom Probekérper aufgenommenen Salzionen darstellt. Die 
jetzt geschilderte neue Erscheinung beweist jedoch, dass sich im Inneren der 
Poren auch noch ein anderer Vorgang abspielt. Dieser andere Vorgang besteht 
darin, dass auch im Inneren der Poren, wenn sich ein grosser Teil der Lésung 
von dort entfernt hat, eine Kristallisation erfolgen kann. 

Dieses Stadium der Kristallisation kann man sich folgendermassen vor- 
stellen: es sei angenommen, dass sich ein so grosser Teil der Lésung entfernt 


Abb. 2. Unterteil des in Abb. 1 gezeigten Fayencetellers mit einer Reihe kleiner Spriinge. 


hat, dass sich an den Innenwdnden samtlicher Poren eine adsorbierte Schicht 
von einer Dicke von nur wenigen Molekilen befindet. Man hat es also in diesem 
Falle — wie bereits oben erwahnt — mit leeren Réhrchen zu tun, in denen 
die Lésung nur in Form einer diinnen, an der Wand adsorbierten Schicht vor- 
handen ist. Das Réhrensystem der Poren stellt ein grosses Labyrinth dar, 
das vielerlei Verbindungen mit der Aussenluft besitzt. Es kann somit Stellen 
geben, wo im Inneren der Poren eine starke Verdunstung einsetzt, weil von 
dem betreffenden Punkt ein kurzer Weg an die Oberflache bzw. ins Freie fiihrt. 
Es kénnen sich aber an einzelnen Stellen im Inneren des Probekérpers auch 
gréssere Hohlriume befinden, von denen ein ganzes Porensystem ausgeht, das 
seine ins Freie fiihrenden Kanile besitzt. An solchen Stellen kann also eine 
intensivere Verdunstung auftreten und als deren Folge das Wachstum einzelner 
Kristalle einsetzen, was schliesslich zum Aufquellen eines Teiles des Probe- 
kérpers filhrt. Die Tatsache, dass es Stellen im Inneren gibt, wo solche grissere 
Hohlraume vorhanden sind, wird durch das System von Aufquellungen ange- 
zeigt, dass dieser Tage an der Unterseite des in Rede stehenden Fayencetellers 
enstanden ist. Die in Abb. 2 sichtbaren Aufquellungen sind regelmdssig am 


—— 
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Rand der Unterseite des Tellers angeordnet, hangen also mit der Herstellung 
des Tellers zusammen. Bei der Herstellung des Tellers diirfte sich namlich ein 
Hohlraumsystem gebildet haben, vermutlich dadurch, dass eine im urspringlichen 
Material befindliche Luftblase zu vielen kleinen Blaschen zerquetscht wurde. 

Die geschilderten Aufquellungen zeigen in augenfilliger Weise, dass sich 
im porésen Stoff, wenn er einmal mit Lésung angefiillt wurde, eine lange Zeit 
hindurch Prozesse abspielen, die auch mit einem bedeutenden Stofftransport 
und einer Verdnderung der Innenstruktur einhergehen. Vielleicht darf man 
hierzu noch bemerken, dass diese Prozesse, die sich an den Porenwinden des 


Abb. 3. Feuerfestes Rohr mit Kristallen und Spuren von Ausbréckelung. 


porésen Stoffes abspielen, niemals zum Stillstand kommen, da an den porésen 
Wanden stets aus der Luft adsorbierbares Wasser zur Verfiigung steht, das 
infolge der Temperatur- und Luftdruckveranderungen — unabhangig von der 
niemals stillstehenden Warmebewegung — standig in Bewegung gehalten wird. 

3. Die hier geschilderte Erscheinung tritt nicht nur bei Fayence auf, 
sondern wahrscheinlich bei jedem porésen Material. Im nachstehenden werden auf 
einigen Lichtbildern Beispiele hierfiir auf verschiedenen porésen Stoffen gezeigt. 

Abb. 3 zeigt ein feuerfestes Rohr, auf dem wollartige NaCl-Kristalle 
sichtbar sind und sich das Material im Stadium eines ziemlich fortgeschrittenen 
Verfalles befindet. Bei Versuchsbeginn bildete die obere Offnung des senkrecht 
aufgestellten Rohres einen regelmassigen Kreis, doch sind jetzt dort infolge 
der Abbréckelung ziemlich grosse Stiicke ausgebrochen. An der Seite des 
Zylinders sind unterhalb der Kristalle auch die starke Schatten werfenden 
Vertiefungen zu sehen. 

In Abb. 4 wird ein gewoéhnlicher Blumentopf gezeigt. Nach nicht langer 
Zeit brachen aus dem Topf grosse Stiicke aus und traten, wie aus Abb. 4/a 
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erstchtlieh betrachtliche Spriinge auf. Neben den auf dem Probekérper befind- 
lichen Kristallen ist ebenfalls die aufgerauhte Oberfliche des Vasenkérpers 


zu sehen. 
In Abb. 5 wird eine kleine emaillierte Blumenvase bzw. eine kleine Schale 


vorgefiihrt. Auch diese Gefisse bestehen aus grobkérnigem Material, wie es von 
den Tépfern auf dem Dorfe zur Herstellung von billigen Gebrauchsgegenstanden 


b) 


Abb. 4, Gewéhnlicher Blumentopf mit Spuren von Aufquellungen und Spriingen. 


benutzt wird. Am Rande der Schale ist gut ersichtlich, wie sich ein grosser 
Teil zusammen mit der Glasur in die Héhe gehoben hat. 

In Abb. 6 ist ein Tiegel aus Aluminiumoxyd zu sehen, der eigens fiir unser 
Institut zum Schmelzen und Ziichten von Alkalihalogenidkristallen angefertigt 
wurde. Indessen konnten die Tiegel nicht fiir diesen Zweck verwendet werden, 
da ihr Porendurchmesser so gross war, dass das geschmolzene NaCl oder andere 
Salze aus ihnen ausflossen. So wurden einige noch unbenutzte Tiegel zur Zichtung 
von Nadelkristallen erprobt. (Wir mussten nimlich systematisch porése Materia- 
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lien suchen und an solchen, von der Industrie erhiltlichen Gegenstinden Proben 
vornehmen.) In Abb. 6/a sind die Salzkristalle sowie die Aufquellungen und 
einige Lécher zu sehen. Abb. 6/b zeigt einen Tiegel, von dem das Salz abge- 
waschen wurde, damit die Verletzungen des Probekérpers gut sichtbar werden. 
Hier ist erkennbar, dass die fehlenden Teile so aussehen, als ob sie von unten 


Abb. 5. Griinglasierte Schale mit einem aufgequollenen Teil. 


her von einem einzigen Punkt aus ausgebrochen worden waren. Abb. 6/c zeigt 
einen kleineren Oberflachenteil mit den Salzkristallen und den Spriingen in einer 
starkeren Vergrésserung. Aus allen diesen Bildern geht hervor, dass auch die 
an der Oberflache enstandenen Salzbliiten einzelne zentrale Rosetten bilden, 
wodurch sie anzeigen, dass an einigen Punkten der Stoffnachschub grdésser ist. 

Abnliche Beobachtungen konnten auch an Bimsstein gemacht werden. 
Auch hier riefen die Salzkristalle Ausbréckelungen hervor. 

4. Die zerstérende und verwitternde Wirkung des gefrierenden Wassers 
in der Natur ist bekannt. Hier handelt es sich also darum, dass diese Erscheinung 
auch bei anderen Stoffen, u. zw. bei der Auskristallisation von gelésten Substan- 
zen auftritt, und es unterliegt keinem Zweifel, dass dies in der Natur haufig 
vorkommt. 

Anlasslich eines Besuches in unserem Institut ausserte Prof. J. D. Bernal 
(London) die Meinung, dass die hier mitgeteilten Beobachtungen eine Erklarung 
fir die allmahliche Verwitterung und Zerbréckelung der Baumaterialien bieten. 
Er fiigte noch hinzu, dass man auf Grund systematischer Untersuchungen 
hoffen kénne, eine Methode zur Verhinderung dieser zerstérender Wirkung 
zu finden, was von grosser praktischer Bedeutung ware. Es scheint, dass man 
nach der Erforschung der Einzelheiten des in den Materialien vorkommenden 
Kristallwachstums auch daran denken kann, entsprechende Methoden eben zur 
Verhinderung der Ausbildung von zerstérenden Kristallen auszuarbeiten. Fiir diese 
Anregung sei auch an dieser Stelle Prof. Bernal der beste Dank ausgesprochen. 
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Abb. 6. Tiegel aus Aluminiumoxyd: a) mit Kristallen und Spriingen, b) dasselbe nach Abwaschen 
der Kristalle, c) ein Teil davon vergréssert. 


Bei den hier beschriebenen Versuchen .waren Adjunkt F. Fathy und 
Laborantin K. Kuch behilflich. Die Kosten der’ Versuche wurden von der 
Ungarischen Akademie der Wissenschaften getragen. 


PASPYLJAIOWJEE TEMCTBYUE PACTYLUUX KPHCTAJIJINTOB 
B TMOPHCTbIX MATEPHAJIAX 
3. bIOSIAM 
Pe3tome 

B NOpucTbIX MaTepHalax, KaK HaNpumep depdop, dasuc u Kepamnyeckue MaTepyhalibl 
nocie BUMTHIBAaHHA pacTrBopa NaCl BpipalMBaiwTcA MrObYaTbIe UM IWePCTAHO-BOOKHMCTHIC 
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MecTax OOpasyl0oTCA pasOyxXaHHA, OCbINAHHA, a TaK)Ke PaspbIBHI. 
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BEITRAG ZUR UNTERSUCHUNG 
DER REKRISTALLISATIONSVORGANGE 


Von 
Z. MORLIN 


INSTITUT FUR EXPERIMENTELLE PHYSIK DER TECHNISCHEN UNIVERSITAT 
FUR DIE BAUINDUSTRIE, BUDAPEST 


(Vorgelegt von Z. Gyulai. — Eingegangen: 25. VI. 1954) 


1. Es wurde der Rekristallisationsmechanismus von NaCl-Pastillen mikroskopisch 
untersucht. 
2. Es bilden sich verhaltnismissig grosse und gut reflektierende Kristalle, deren Flachen 
eine feinere Struktur aufweisen. 
3. Versuche zur Deutung der beobachteten Feinstruktur wurden unternommen. 
4. Die Frage der Korngrenzen wurde kurz besprochen. 


In einer friheren Arbeit wurde tber Rekristallisationsversuche an NaCl- 
Pastillen, die aus feinem NaCl-Pulver unter hohem Druck hergestellt wurden, 
berichtet [1]. *Es konnten Erholungserscheinungen auf elektrischem und 
optischem Wege ausgemessen werden, ferner wurden mikroskopische Korn- 
zahlungen vorgenommen. Es zeigte sich, dass das durchschnittliche Mass der 
Rekristallisation vom angewandten Pressdruck und von der Presstemperatur 
abhangt. Es wurde festgestellt, dass schon eine kurze Pressdauer (die Herstellung 
der Pastillen dauerte 60 Sekunden) zu einem bedeutenden Kristallwachstum 
fuhrt. Héhere Presstemperatur begiinstigt den Vorgang. Diese Untersuchungen 
wurden fortgesetzt. Vom Standpunkt des Kristallwachstummechanismus erhebt 
sich nun die Frage, auf welche Weise die vielen verschieden orientierten einzelnen 
Kristalle zu einem einzigen grésseren zusammenwachsen. 

1. Wenn man nun das Kristallwachstum aus der festen Phase untersucht, 
so ist die Tatsache zu beriicksichtigen, dass sich bei der Rekristallisation die 
Wachstumsumstande anders gestalten als beim Wachstum aus der Lésung, der 
Schmelze oder der Dampfphase. Im letzteren Falle wird der wachsende Kristall 
von freien, sehr beweglichen Ionen oder Ionenpaaren umgeben, im Falle der 
Rekristallisation geht man hingegen von einem Kristallhaufen aus, in dem die 
Ionendichte annahernd dieselbe ist wie beim Einkristall. Der schnelle Verlauf 
des Kristallwachstums bei den Pastillen ist nun dadurch erklarbar, dass die 
Ionen hier in entsprechender Dichte vorhanden sind, und die einzelnen Ionen 
sich nur einige Ionenentfernungen zu verriicken brauchen, damit die notwendige 
geometrische Ordnung hergestellt wird. 

Das Ausgansgmaterial wurde ebenso behandelt, wie frither. Das im Achat- 
morser zerriebene NaCl-Pulver wurde in einem Druckbereich von 1000—20 000 
kg/cm? mit einer Presstemperatur von 20—750 °C zu Pastillen gepresst. Die 
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Abb. 1, Flachen rekristallisierter Einkristalle an frisch gespalteten Pastillenoberflachen. 
la. Pastille Nr. 26. Pressdruck 3000 kg/cm?, Presstemperatur 750° C 200x. 
2b. Pastille Nr. 27. Pressdruck 5000 kg/cm?, Presstemperatur 750° € 200x. 
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fertigen Pastillen wurden auf dreierlei Weise untersucht. Auf den Pastillen- 
oberflachen wurden gut reflektierende Flachen beobachtet. Wenn diese Flichen 
mit Oberflaichen frisch gespalteter Einkristalle verglichen werden, sieht man, 
dass die gut reflektierenden Flachen rekristallisierte. Einkristalle in der Pastille 
sind [2]. Dafiir sprechen jene feinen Linien, die nur bei Einkristallen beob- 
achtet werden kénnen. Abb. la und lb zeigen solche Flachen, die bei jedem 
Druck, und bei jeder Presstemperatur zu beobachten sind. Bei Anwendung 
héheren Druckes und héherer Presstemperatur nimmt jedoch die Zahl und der 


Abb. 2. Pastille Nr. 20. Flache eines rekristallisierten Minicrievabes an einer frisch gespalteten 
Pastillenoberflache. Pressdruck 3000 kg/cm?, Temperatur 300° C 800x. 


Umfang dieser Flachen zu, dabei werden die rekristallisierten Kérner nicht nur 
grésser, sondern es bilden sich auch bessere Kristalle, hierbei spielt die Press- 
temperatur eine wichtige Rolle. Abb. 2 zeigt die Flache eines rekristallisierten 
Kristalls aus einer Pastille, die mit 300° C und 3000 kg/cm? Druck hergestellt 
wurde. Die Flache ist nicht so gut wie die Flache auf Abb. 1a oder 1b. Hingegen 
zeigt Abb. 3 einen grossen Kristall, der bei einer Presstemperatur von 750° C 
mit 6000 kg/em? Druck hergestellt wurde. 

Zur Untersuchung der Korngrenzen wurden die Pastillen nach ooo 
Spalten poliert, dann in 94%-en Alkohol gedtzt. Abb. 4 (4a und 4b) zeigt zwei 
Pastillenflachen, die poliert und dann geatzt wurden. Diese Behandlung brachte 
die Korngrenzen, d.h. die Konturen der rekristallisierten Einkristalle gut 
zum Vorschein. Es ist bermerkenswert, dass die einzelnen Kristallite der Pastille, 
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die mit 3000 kg/cm? und 20° C gepresst wurden viel kleiner sind (Abb. 4a) 
als die Kristalle jener Pastille (4b), die mit 10 000 kg/cm? und 100° C her- 
gestellt wurde. 

2. Zu den weiteren Untersuchungen wurden dieSpaltflachen der Pastillen 
gleich nach dem Spalten in 94%-en Alkohol geatzt. Nach einer solchen Behand- 
lung zeigen die ursprimglich strukturlosen Kristallflachen in der Pastille eine 
feine Struktur (Abb. 5, 6, 7, 8). Zur Herstellung der diesbeziiglichen Mikro- 
photographien ist zu erwahnen, dass die Spaltflachen der Proben nur selten 


Abb. 3. Pastille Nr. 29. Rekristallisierter Kristallit nach einstiindigem Pressen. Der Kristall 
wurde aus der Pastille herausgebrochen. Pressdruck 6000 kg/cm?, Presstemperatur 750° C 65x. 


ganz eben waren, daher konnte immer nur ein Teil der beobachteten Pastillen- 
oberflachen scharf eingestellt werden. Die abgebildeten Flachen sahen urspriing- 
lich vor der Atzbehandlung ganz so aus wie jene, die in den Abb. la und Ib 
wiedergegeben sind. Dass es sich hierbei um Einkristalle handelt, ist durch das ein- 
heitliche Reflektieren der Flachen erkennbar. Nach der Atzbehandlung zeigen 
diese Flichen feine Atzfiguren, welche die feinen Verschiedenartigkeiten in der 
Struktur der Oberflichen zum Vorschein bringen. Die Atzfiguren bestehen aus 
kleineren und grésseren quadratischen Flachen, deren Diagonallinien scharf 
hervortreten. Die quadratischen Flachen sind zueinander parallel gelagert ; 
also sowohl die Seitenlingen wie die Diagonallinien sind zueinander parallel. 
Die ganze Flache eines Kinkristalls ist mit diesen einheitlich gelagerten quadra- 


> 
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Abb. 4. Polierte und geadtzte Pastillenoberflachen. 
4a. Pastille Nr. 23. Pressdruck 3000 kg/cm?, Presstemperatur 20° C 200x. 
Ab. Pastille Nr. 24. Pressdruck 10,000 kg/cm?, Presstemperatur 100° C 200x. 
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tischén’ Figuren bedeckt. In den photographischen Aufnahmen (z. B. in Abb. 5a) 
sieht man neben einigen scharf hervortretenden Quadraten teilweise ver- 
schwommene Linien, die wegen Unebenheiten der Flache nicht scharf hervor- 
treten, und die ahnen lassen, dass die ganze F'lache diese feine Struktur enthalt. 


Abb. 5a. Pastille Nr. 25. Oberfliche eines neu gebildeten Kristalls nach Atzung der Pastillen- 
oberflache. Die Feinstruktur wird sichtbar, es zeigen sich Wiirfelflachen mit Querstreifen 
und Konturen eines Netzssystems. Pressdruck 3000 kg/em?, Presstemperatur 700° C 480x. 


Abb. 5b. Zeichnerische Abbildung der Aufnahme 5a zur richtigen Orientierung. 
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In Abb. 5b, die die in Abb. 5a wiedergegebene mikroskopische Beobachtung 
” zeichnerisch darstellt, wurden diese nicht allzu dominierenden Linien hervor- 
gehoben. Eine ahnliche Struktur zeigt Abb. 6a und die dazu gehérende Zeichnung 
Abb. 6b. Abb. 7 und 8 zeigen Einkristalle, die dieselbe Struktur aufweisen. 


Abb. 6a. Pastille Nr. 42. Weiteres Beispiel der Feinstruktur der neugebildeten Kristallite. Das 

Bild zeigt einen Flachenteil eines rekristallisierten Einkristalls in der Pastille. Neben den scharf 

eingestellten Konturen sieht man weitere Konturen ahnlicher Form. Pressdruck 5000 kg/cm?, 
Presstemperatur 750° C 780x. 


i 


Abb. 6b. Zeichnerische Abbildung der Aufnahme 6a zur richtigen Orientierung. 
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Die besprochenen Abbildungen sind ahnlich einer von Gyulai beobachteten 


Wachstumsfigur [3], die in Abb. 9 wieder abgedruckt ist. Diese Wachstums- 
figur stellt die Ausbildung eines Kristallskeletts dar, das aus einer Lésung bei 
sehr raschem Wachsen enstanden ist. Das Wachsen begann am Rande des 
Tropfens. Man sieht, dass einzelne Kristalle sich zueinander ordnen, dabei 
wird die diagonale Richtung der quadratischen Kristalle beibehalten. Man 
sieht ferner, dass jedes Eck ein wachstumsfahiger Punkt ist an den neue Kristal- 
lite mit ihren Diagonalrichtungen anwachsen. 

Die Ahnlichkeit der Abb. 9 mit den vier fritheren Abbildungen legt den 
Gedanken nahe, dass wir in den Abb. 5, 6, 7 und 8 die Wachstumsfiguren der 


Abb. 7a. Pastille Nr. 40. Abbildung einer typischen Wiirfelflache mit Querstreifen. Die beige- 
fiigte Zeichnung zeigt als Fortsetzung in der einen Diagonalrichtung eine ahnliche Flache. An 
der Photographie ist nur die linke obere Spitze der letzteren sichtbar. Die ganze Flache bildet 
ein Netzsystem, deren einzelne Elemente sich liickenlos zuasammenfiigen. Pressdruck 3000 kg/cm2, 
Presstemperatur 750° C 780x. 
Abb. 7b. Zeichnerische Abbildung der Aufnahme 7a. 


rekristallisierten Einkristallchen vor uns haben. Selbstredend kénnen bei der 
Rekristallisation nicht einzelne Kristallreihen enststehen, da hier der ganze 
Raum voll mit Ionen ist, also an jedes Korn von allen Seiten andere Kristallite 
anwachsen. Man kénnte sagen, dass der in Abb. 9 dargestellte Vorgang, der 
sich quasi immer in linearer Richtung abspielt, bei der Rekristallisation zwei- 
dimensional sichtbar wird. In Wirklichkeit spielt sich die Rekristallisation 
dreidimensional ab, davon ist aber in der Spaltflache nur der zweidimensionale 
Teil sichtbar. Bei eingehender Untersuchung der photographischen Aufnahmen 
sieht man, wie dies schon hervorgehoben wurde, dass die ganze Flache eines 
Einkristalles mit den Spuren quadratischer Figuren iibersat ist, die aber in, den 
Aufnahmen nicht scharf genug hervortreten. Schliesslich ist es nur eine Frage 
der Beobachtungstechnik, die besprochene Struktur in jedem einzelnen Ein- 
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kristall einer rekristallisierten Pastille nachzuweisen. Man kann folglich mit 
einer gewissen Anniherung behaupten, dass das Wachstum eines Einkristalles 
in einer Pastille von einem Kern ausgehend schrittweise vor sich geht. Die Re- 
kristallisation geht vom ersten Kern aus, im weiteren wachsen dann an jedes 


Abb. 8a. Pastille Nr. 38. Oberflache eines neu gebildeten Kristalls nach Atzung der Pastillen- 

oberflache. Die Feinstruktur zeigt Hauptwachstumsrichtungen, Fortsetzung einer Diagonale, 

Flachen mit Querstreifen und Zeichen schichtenweisen Wachstums. Pressdruck 5000 kg/cm?, 
Presstemperatur 750° C 780x. 


Abb. &b. Zeichnerische Darstellung der Aufnahme 8a zur richtigen Orientierung. 


Kristallchen gleichgerichtete Nachbarkristallite liickenlos und gleichorientiert 
an. Der Vorgang spielt sich auch hier sehr rasch ab. Der Pressdruck wurde nur 


2 Acta Physica IV,3. 
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einige Minuten angelegt, es wurden aber auch Proben mit einer Pressdauer 
von einigen Sekunden hergestellt. Die Pastillen sehen in jedem Fall gleich aus. 
Die fiir spiter geplanten quantitativen Versuche diirften in dieser Richtung 
vermutlich unsere Kentnisse vertiefen ; soviel kann jedoch schon jetzt behauptet 
werden, dass die Ausbildung eines Einkristalles in einer Pastille in raschen, 
aufeinander folgenden Schritten vor sich geht, und zwar so, dass die kleinen 
quadratischen Gebilde (wobei es sich um dreidimensionale Gebilde handelt) 
sich sprungartig nebeneinander reihen. Dieses Bild entspricht — im drei dimen- 
sionalen Vorgang — der Kossel—Stranskischen molekulartheoretischen Be- 
trachtung. Beobachten wir jene Wachstumsfiguren von Gyulai, die Kristallisa- 
tion aus Dampf und Lésung zeigen [3] und vergleichen wir dieselben mit den 
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Abb. 9. Diese Abbildung ist der Arbeit Z. Gyulais, Z. Kristallogr. (A) 91/ 142, 1935, 
Seite 146 entnommen. 


Abb. 5, 6, 7 und 8, die die Wachstumsfiguren bei der Rekristallisation darstellen, 
so ist die Ahnlichkeit der Abbildung auffallend. Man kann diese Ahnlichkeit 
dahin deuten, dass das Wachstum in jeder Form (Lésung, Dampf, fester Zustand) 
nach denselben Regeln vor sich geht, und zwar gemiiss den Kossel—Stranski- 
‘schen Regeln. Wir sind uns dabei der Tatsache bewusst, dass sich die Kossel— 
Stranskische Theorie auf einzelne molekulare Vorginge bezieht, wahrend das 
Wachsen bei der Rekristallisation schichtenweise vor sich geht. Jedoch muss 
hier in Betracht gezogen werden, dass, wie Gyulai gezeigt hat, die Kosselsche 
Regel auch das schichtenweise Angliedern der Molekiile wiedergibt [4]. 

3. Es sei noch eine Beobachtung kurz besprochen. In den Abb. 4a und 4b 
sind Pastillenflachen wiedergegeben, an denen die einzelnen Einkristalle, in der 
schon besprochenen Weise sichtbar gemacht wurden. Zwischen den einzelnen 
Kristalliten sieht man kleine Zwischenriume. Man sieht, dass die eingeinen 
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Kristalle nicht unmittelbar aneinander grenzen, sondern dass zwischen den 
einzelnen Kristalliten ein gewisser Abstand vorhanden ist. Dies ist bei jedem 
Druck und bei jeder Temperatur zu beobachten. In der diesbeziiglichen Literatur 
wird nun die Frage vielfach besprochen, dass die Zwischenraume eines poly- 
kristallinen Aggregates aus amorphem Material bestehen, bzw. dass den Raum 
zwischen den einzelnen Kérnern unter Umstinden eine etwa 100 Atome um- 
fassende Schicht ausfiillt, die wegen der in ihr herrschenden starken Spannung 
als quasi amorph betrachtet werden kann [5]. In den Abb. 4a und 4b sieht 
man die einzelnen Einkristalle je einer Pastille. Zwischen den einzelnen Kérnern 
sind die erwihnten Zwischenraume sichtbar. Zur Entstehung der Zwischen- 
raume sei nun folgende Bemerkung gemacht. Aus eigenen noch nicht veréffent- 
lichten Dichtemessungen wissen wir, dass die Dichte der Pastillen geringer ist 
~als die der Einkristalle. Die Dichte der mit grésserem Pressdruck hergestellten 
Pastillen ist grésser als die Dichte jener Proben, die mit geringerem Druck 
hergestellt wurden. Bei dem Vorgang der Rekristallisation, also im Bereich 
der einzelnen Einkristalle, dichtet sich das Material zusammen, demzufolge 
zwischen den einzelnen Kérnern kleine Zwischenraume entstehen und zwar 
dort, wo infolge der Schrumpfung des Materials die innere Spannung gleich 
Null wird. Dort, wo die Spannung aufhért bzw. unter eine gewisse Grenze gefallen 
ist, hért die intensive innere Beriithrung zwischen Einkristall und Umgebung 
auf, und das Wachsen des Einkristalls kommt zum Stillstand. In den Zwischen- 
raumen bleibt also eine grosse Anzahl] Kristallsplitter angehauft, die mitein- 
ander nur locker verbunden sind und den Raum lose ausfillen. Da nun das 
Atzen zunichst die schwiachsten Teile der Pastillenoberflachen angreift ; werden 
zuerst die lockeren Splitter aufgezehrt, die zwischen den rekristallisierten 
K6rnern eine Schicht bilden, welche zu keinem der rekristallisierten Kristalliten 
gehért. Daher die beobachteten Zwischenraume. Aus diesen Betrachtungen 
folgt, dass das Vorhandensein einer amorphen oder qusai-amorphen Schicht 
zwischen den Korngrenzen bei NaCl-Pastillen als unwahrscheinlich betrachtet 
werden kann. 

Die Arbeit wurde im Institut fiir Experimentelle Physik der Technischen 
Universitat in Budapest mit Mitteln durchgefirt, die die Ungarische Akademie 
der Wissenschaften zur Verfiigung stellte. Aufrichtiger Dank sei dem Leiter 
des Instituts, Herrn Prof. Gyulai fiir die Anregung der Arbeit und stetige 
Leitung derselben ausgesprochen. Ebenso danke ich Herrn Mechaniker Karl 
Wallentin fir die Herstellung der Pressformen. 
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BUDAPEST 


(Vorgelegt von K. F. Novobdtzky. — Eingegangen: 9. VII. 1954.) 


Durch Kombination verschiedener Methoden wird das elektromagnetische Feld mit Hilfe 


_ einer einzigen komplexen Grosse beschrieben und gequantelt. Die mit der Erfahrung iiberein- 


stimmenden Resultate lassen sich anschaulich leicht deuten. 


Es ist bekannt, dass die Quantelung des elektromagnetischen Feldes 
mehrere Schwierigkeiten in sich birgt. Zu deren Eliminierung wurden ver- 
schiedene Versuche gemacht. (S. z. B. Heisenberg—Pauli [1], Fermi [2], 
Chang [3] u. s. w. ) 1938 gab Novobdtzky [4] eine sehr erfolgreiche Methode an, 
die im wesentlichen darauf beruht, dass das skalare Potential im Fourier-Raum 
durch das Vektorpotential eliminiert wird. Vor kurzem gelang es Green und 
Wolf [5] mittels komplexer Kombinationen des Vektorpotentials, die Lagrange- 
Funktion des elektromagnetischen Feldes auf eine dusserst einfache Form zu 
bringen. Dieses Verfahren ist aber nur im Falle eines reinen Strahlungsfeldes ' 
anwendbar. Im nachstehenden werden wir zeigen, dass das Green und Wolfsche 
Verfahren mit der Methode von Novobatzky kombiniert es erméglicht, den 
kanonischen Formalismus auch im Falle des allgemeinen, auch Ladungen 
enthaltenden elektromagnetischen Feldes in einer erstaunlich einfachen Form 
auszuarbeiten. Die nachfolgende Quantelung bietet keine Schwierigkeiten. 

Ohne Beniitzung der Lorentzschen Bedingung kann man die Maxwellschen 


Gleichungen in folgende Form bringen : 


AD + tae Pe Y= — 4x0, (1) 
c Ot 
2 1 : An , 
AYA eee, grad (div % + — P| = hee; (2) 
ea. i é c c 


Entwickeln wir die Potentiale und die Ladungs- und Stromdichte in 


ablicher Weise in eine Fourier-Reihe : 


D = Qf elt (3) 
k 
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Ses Qn XY aeit® (4) 
is 
OS Qa get (5) 
k 
{= 21? ize ite, (6) 
k 


Fir das Folgende ist eine spezielle Komponentendarstellung des Koeffi- 
zienten 4, sehr zweckmassig. Wahlen wir zwei aufeinander und auf £ ortho- 
gonale Einheitsvektoren so, dass ¢,*, ¢,* und f ein Rechtsystem bilden. In 
diesem System hat 4, die folgende Gestalt : 


f 
irs OO) tae Cora ra (7) 


Unsere erste Aufgabe besteht in einer Umformung des skalaren Potentials 
im Fourier-Raum. Das kénnen wir nach Novobatzky erreichen. Aus (1) haben 
wir namlich 


also statt (3) schreiben wir : 
Ba O-1u2 yh 45) 4 Ae | (3’) 
k } ck2 . k2 


Als nichstes folgt die Eliminierung von a,", des longitudinalen Anteils 
in 4. Dazu bietet die Eichinvarianz der Feldgleichungen eine Méglichkeit.. 
Diese sind namlich invariant gegeniiber der Transformation 


P+ 6446 
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: U—> W— grad G, 


G = QD Y gy elt 
k 


sei. Deshalb haben wir mit (7) 


A S| a ef + ob eh + (af — thay) Le 
k ¥ 


Ks sei 


1 
k 
Sk S= = 743; 
tk 
/ 
dann treten in der Reihe von nur die transversalen Komponenten auf, nam- 


lich es wird jetzt 


ay, = at ef + ak ef. (8) 


Wir bemerken, dass die obige Transformation von 2 in (3’) auch schon 
die richtige Transformation von ® zur Folge hat, wie man leich nachpriift. 

Nach Eliminieren des longitudinalen Feldes bekommt ® selbst eine ein- 
fache Gestalt. Nach (8) ist 


(fax) = 0, (9) 


und deswegen kann man das skalare Potential durch die gegebenen Quellen 
ausdriicken : 


An 0 
o-girny™, =(2av. 3””) 
py ae ( 


Nun haben wir also zwei unabhangige Veranderliche a, und a,", die sich 
wie Vektorkomponenten transformieren. Wegen der Realitét von Y gilt 


a; = Cy, (10) 
das heisst 
—k = a‘* 
ay fit (11) 
ag — a9 5 


falls wir verabreden, dass wir zu —f£ ein Linkssystem zuordnen, also, dass 
a ke. Sea Si 
Gi eis ex es 


ist. 


_ Stromdichte : 


“s 


=~ Orticteis (11) und (12) lauten die inversen Formeln : 


» 


Se ce be Ne A 
at = gre + a1), 


Wir konstruieren nun aus den a,-s eine einzige Feldfunktion. Es sei. 


p= 22S ay ote, C4) 


Diese Funktion charakterisiert das Feld eindeutig, da aus ihr riickwarts die 
elektromagnetischen Potentiale rekonstruierbar sind. 
Bevor wir weitergehen, zerlegen wir auch die Fourier-Komponenten der 


ix = ji ef +73 es +75 e5. > (15) 
Es sei : 
Ay = ji + 8 — (16) 


Mit Hilfe der Realitat von j haben wir in analoger Weise 


repel ory Pepa | 
Ji 9 (Ak + ake | 


(17) 
1 : 
ja — — (Ay, An x) 
2 
Definieren wir noch die Grésse 
T= 2-1? SA, eft (18) 
fir 


Nach diesen Vorbereitungen kénnen wir die vollstandige Lagrange- 
Funktion konstruieren. Bekanntlich ergeben sich die Feldgleichungen (1), (2) aus 


Les OF r i 
= | a + q Merad D+ grad? ® — ror? a | — {A — o® jay. 
2 


(19) 
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Hierin die Reihenentwicklungen (3”), (4), (5), (6) eingesetzt, erhalten 
wir mit Riicksicht auf die Realitatseigenschaften der F ourier-K ompunenten 
sowie mit Hilfe von (8) bzw. (15) und (13) bzw. (17), ferner den Definitionen 
(14) bzw. (18) das dusserst einfache Ergebnis 


iF “Al 


— pp* — grad y grad y* | + =~ (le" + I*p) — 


8x | c? 


1 E 
—_O- Uy: ans Se es (20) 


Diese Lagrange-Funktion stimmt, wenigstens dusserlich, mit der Lagrange- 
Funktidn des komplexen skalaren Feldes mit der Ruhmasse Null iiberein, abge- 


sehen von der letzten unwesentlichen additiven Konstante. Die Feldgleichungen 


lauten : 
c 
und (21) 
pt gy* — An Ts, 
c 
Die kanonisch konjugierten Feldgréssen sind : 
T= 22 os 
dp 82c? 
(22) 
PTR OL ee | 
aqg* Bac? 
Die Hamilton-Funktion lautet also: 
] Sk Sp 
Be | Nie (a =p *p+grad p grad g* | aay (Ig* + Ip) av a vee k Sk 2 (23) 


Das letzte, von den Feldgréssen unabhingige Glied ergibt die Coulombenergie. 


Es ist namlich 


sue Q-2 0 (r) e— itt dV, 
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und daher 


sat : eif(e’—r) ; | 
2a >. orn =s 2x2 | loft) o(t’) Dai are a dVdV’ = 


k 


= on (= dfo(t) o(t’) dVdV" = Sarg dVadv’, 


w. Z. b. w. 

In ganz ahnlicher Weise kénnen wir auch beliebige andere integrale 
Feldfunktionen, z. B. den Feldimpuls, mit Hilfe von q allein eindeutig kon- 
struieren. Als sehr interessant erweist sich z. B. die Berechnung des Spins. 
Nach der allgemeinen Definition (s. z. B. [6]) ist 


S2= \a ie ydV, (24) 


worin [12 den Infinitesimaloperator der Drehung um die Z-Achse bedeutet. 
Zunichst wollen wir diesen berechnen. Da aj und aj Vektorkomponenten 
sind, transformieren sie sich bei der genannten Drehung wie folgt : 


aj = ai + agde, | 


oF ta wees 
Also ist 
a, = at + iak + (ak —iak)de = a, — ia;,de, (26a) 
und ahnlich 
a, =a, +-ia; ds. (26b) 


Wenn wir jetzt im f-Raum bei jeder einzelnen Partialwelle eine Drehung um 
die f-Richtung um ein und demselben Winkel d¢ vollfiihren, so transformiert 
sich selbst folgendermassen : 


y’ = yp — ioe, | 
(27) 
gt’ =g*+ide . 


und 


Ubrigens lauten diese Transformationsgleichungen in integrierter Form : 


yp’ =e", 
(27) 


g*’ = pte'® P 


ae 
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Konstruieren wir aus p und y* eine zweireihige Matrix 


ha 

oo | yr |” 

so sehen wir, dass 
: —10 
v=ptil ay | wae, 

also 

2 |—10 
ist. Nach (24) ist also: 

4 
Sk =| IT IT* | pe oe dV = i| o* — II) dV. (29) 


Wir bemerken, dass die X und Y Komponente des Spins Null ist. Wenn wir 
namlich die Drehung nicht um die f, sondern’z. B. um die es Achse vollziehen, 
so ist 


Daraus folgt aber, dass diejenigen Matrixelemente von [*°, die auf m und q* . 
wirken, gleich Null sind. 

Diese Berechnung des Spins steht in Einklang mit den Resultaten von 
Pauli [7] und Shirokov [8], denen zufolge der Spin eines Teilchens mit ver- 
schwindender Ruhmasse immer in die Fortpflanzungsrichtung hinweist. 

Da das g-Feld komplex und gegeniiber Eichtransformation erster Art 
invariant ist, erhebt sich die Frage, ob eine «Feldladung» existiert. In bekannter 
Weise ist 


e=— in| (¢* p—or* dV, (30) 


dies stimmt aber, abgesehen von einem Proportionalitatsfaktor mit dem Spin 
iiberein, also bekommen wir keine neue Feldgrésse. Dies miissen wir ja auch 
von vornherein erwarten da, wie z. B. aus (27’) ersichtlich ist, bei uns die 
Drehungsformation mit der Eichtransformation tbereinstimmt. 

Es ist noch bemerkenswert, dass im vorliegenden Fall der Spin selbst 
einem Erhaltungsgesetz gehorcht, wenigstens in Abwesenheit von Quellen. 
Das hangt damit zusammen, dass die. Feldgleichungen gegeniiber der Trans- 


Y 
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formation (27) auch dann invariant sind, wenn wir die Argumente der Feld- 
funktionen nicht transformieren. 
Im nachfolgenden wollen wir noch kurz die Quantelung des Feldes erledigen. 
Die Vertauschungsrelationen lauten : 


UT (0), oe) = 3 Le 9] = 2 e—8). (31) 


Diese sind erfiullt, falls 
fay, | = a Ox Bate, 
L 


Es sei 


= Ay. (32) 
sodann ist 


h 
[vies One] = 2 xn. (33) 


U 


Bestimmen wir jetzt z. B. die Eigenwerte des Strahlungsanteils von H. 
Es wird nach (23), (14) und (32) 


* k2 * 
Hy, = 3} \(82c®) ye ye + oe Ay Ay, | 
a 


: ie 


Zur Darstellung der Operatoren a, und y, fiihren wir die Operatoren c, und d, 
ein, die den Vertauschungsrelationen 


es eure. (di, dx] = Sux (34) 


geniigen (die gemischten Kommutatoren verschwinden). Mit diesen ist 


(35) 


und daher 
hikle * » * * 
H, ee eee ea ) (¢, Ce + dj, dx + cx Cy + dy d;). 


+ k 


ee 
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Nach (34) ist aber 


. * x vf 
der Eigenwert von c; c, gleich Ni. 
* 


« ae <a UL. Cae « Nae, 
« « Rud aK NY, 
« « @ + cd, dj, « Neer 
und deshalb 
Ho = ) bl k| (NY + NP) + Bp: (36) 
ke 


ona a (2) 9. : : : 
Hier sind Nj.) und N;,”? die Besetzungszahlen von zwei verschiedenartigen Photo- 
nenzustanden. 


Jetzt folgt die Quantelung des Spins. Gemiass (29), (14), (22) und (35) 


haben wir 


ee 2G ck — di dx + cect — dy d}), 
k 


also ist 
S; =D) 5 (Ni — NP). (37) 


Wir sehen, dass der Spin geradzahlig ist, wie es sein muss. Ferner bemerken 
wir, dass die obengenannten zwei verschiedenartigen Photonensorten sich eben 
im Vorzeichen des Spins unterscheiden. Weiter kénnen wir noch zeigen, dass 
diese «Polarisation» der Photonenzustande gerade mit der optischen Polari- 
sation der Partialwellen iibereinstimmt. Es sei namlich eine links zirkular- 
polarisierte Welle durch 


" var 
a; = u, cos | k | ct, 


k 
2 


. | | 
az = u, sin | k | ct 


dargestellt. Dann ist nach (12) b. z. w. (32), sowie (35) 


G0 


‘i= / kL Ue Pe 20, 
| Acthc* 


Im Falle einer rechts-zirkularpolarisierten Welle gilt alles mutatis mutandis. 
Daraus sehen wir, dass No die Zahl der rechts-, N,.? die Zahl der links-zirkular- 


polarisietren Zustande bedeutet w. z. b. w. 
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Endlich kénnen wir noch den Feldoperatoren y und g* eine sehr anschau- 
liche Bedeutung geben. Man iberzeugt sich namlich mit Hilfe von (29) und 
der Vertauschungsrelation (31) leicht davon, dass, wenn x einen Zustand 
mit dem Spineigenwert S; bedeutet, so gehért zum Zustand (py) ein um 4 
verminderter Spineigenwert. Ahnlicherweise gilt fir y* das entgegengesetzte. 
Also vermehrt der Operator g* den Spinwert um 1, und ¢ erniedrigt ihn um 
eins. Dies ist schon auch aus der Tatsache zu erwarten, dass, wie gesagt, bei 
uns der Spin- und der «Ladungsoperator» zusammenfallen.* _ 

Zusammenfassend wollen wir herausheben, dass das obige Verfahren zur 
Darstellung und Quantelung des elektromagnetischen Feldes eine bewunderns- 
wert einfache und ibersichtliche Methode ist, mit sehr iibersichtlichen und 
anschaulich leicht deutbaren Resulateten, undzwar ohne Heranziehung von 
lastigen Nebenbedingungen. Natiirlich erhebt unsere Methode keinen Anspruch 
auf Lorentzinvarianz. 

Abschliessend spreche ich Herrn Dozenten G. Marx meinen innigsten 
Dank aus fiir die Anregung vorstehender Arbeit und fiir sehr wertvolle Dis- 


kussionen. ‘ 
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Be: 
2 
: 
; 


? 


CLASSICAL MOTION OF THE NUCLEONS IN A SCALAR 
MESON FIELD 


By ‘ 
G. SZAMOSI and G. MARX 


INSTITUTE FOR PHYSICS OF THE ROLAND EOTVOS UNIVERSITY, BUDAPEST 


(Presented by K. F. Novobatzky. — Received 29. VII. 1954) 


The present paper deals with the classical motion of a point nucleon in a scalar meson 
field. The motion in a static field can be subjected to exact treatment. The calculations corfirm 
Werle’s results obtained for small energy approximation according to which a purely attrac- 
tive (e. g. Yukawa) potential will lead to repulsion if the absolute value of the potential energy 
becomes greater than the proper energy. The existence of a repulsive sphere of a radius app- 
roximately equal to the nucleon’s Compton-wavelength around the nucleon is thereby proved 
without recourse to any special assumption. A study of the course of motion in time leads, 
among others, to the result that the particle crosses the boundary of the repulsive region with the 
velocity of light. In the case of the interaction of several nucleons effects analogous to 
three-body forces take place. 


1. Introduction 


It is a generally accepted view that relativistic effects play no significant 
role under the conditions prevailing in the nucleus. This is based on the fact 
that the binding energy of a nucleon is relatively small compared with the 
proper energy of it. However, no specified statements on this question are 
justified until after an examination of the exact relativistic equation of motion 
of the nucleons in the nuclear field (meson field). At first, the solution must 
be obviously sought for within the framework of classical theory. 

Several authors dealt with the setting up of the equations of motion [1]; 
in addition to the force exerted on the nucleon by the incoming field they deter- 
mined the terms arising from the reaction of self-field. Accordingly, the classical 
relativistic equation of motion is exactly known. Less attention has, however, 
been paid to the solution of the equations, i. e. to the actual course of the space- 
time motion. In discussing the relativistic Hamilton-Jacobi equation Werle [2 | 
recently pointed out that as long as the sum of the kinetic and potential energies 
of a particle is small relative to the proper energy, the motion of the particle 
corresponds to a motion in a potential field with a strong short-range repulsive 
potential added to the common attractive nuclear potential. This result is 
most interesting because if it is confirmed the attractive nuclear force must 
turn into repulsion at a very small distance merely for relativistic reasons. 
This means that the nucleons are surrounded by a small repulsive sphere which 
prevents them from approaching each other unlimitedly. It is known that the 
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existence of such a repulsive sphere was also made probable by Lévy’s quan- 
tum-theoretical considerations [3]. The existence of this short-range repulsion 
seems to be in agreement with experiment, for by the phenomenological 
introduction of it Jastrow and others [4] were able to explain the results of 
nucleon-nucleon scattering experiments and other nuclear-physical observa- 
tions. 

By discussing rigorously the solutions of the equation of motion of 
a nucleon in the classical scalar meson field, this study aims at ascertaining 
whether in a certain region the meson potential gives rise to repulsion for 
relativistic reasons. The problems of nuclear physics are typical many-body 
problems, wherefore the interaction of nucleons too should be investigated 
as a many-body problem. But as no satisfactory solution of the classical 
relativistic many-body problem is known, our examination shall also refer 
to a single point nucleon in a given external meson field. This means that the 
reaction of the self-field will not be taken into account; a procedure, which, 
from the point of view of the problems considered. seems permissible. 


2. The equations of motion with three-dimensional notation 


The equation of the meson field generated by nucleons of nucleon charge g, 
and described by density-function @ is 


19 — #p = dog . (1) 
As is known, the energy-momentum-tensor of the meson field is : 


| eo Me) ay dix ( Op Op 
4a Ox; Ox, 82 | Ox, Ox, 


Lips 


The density of the ponderomotive force acting on the nucleon is 


OTK ) 
= "aS See Ss 


Ox k Ox; 


is 


whence we obtain for the four-vector of the force acting on the point nucleon 


Taking this into account the equation of motion of the nucleon is : 


d Op 
— Mu, = —_—. ; =—=—_—_— 
2 ( ) & ae ° uj uj =—c?, (2) 


a + 
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where Mis the proper mass, u; = dx; / dt the four-velocity of the nucleon, = tel, “ 
c denotes the velocity of light, and t a scalar parameter, which is mostly known 
as «proper time». Introducing the (scalar) potential energy gy = V, we may 
write: 


d OV 
-— (Mu;) = ——_., 3 
dt ( i Ox; ( ) 


Novobdatzky [5] has shown that in the case of a force described by a scalar 
potential energy the rest mass cannot be constant during the motion, for by 
differentiating and thereafter multiplying the equation of motion (3) by ui / c?, 
we obtain : 


V 
Sa M=m+-. 4 
dt c2 dt es (4) 


This permits substituting into the equation of motion the constant nucleon 
mass m known from mass spectrographic measurements : 


[let]. 


Let us rewrite this equation in three dimensional notation. The first three 


components (i = 1, 2, 3) are: 


m+— pa Se 
d 2 = 
Sa) | = “all 1 _ 18 grad V. (6) 
c 


The fourth component (i = 4) multiplied by (—ic) : 


d iad i pro ae . 


#( = ae 
aie. 
Cc 


Let us assume for a moment that V is not explicitely dependent on time (the 


field is static). Then the quantity 


mc2 + V 
yi= He : 


3 Acta Physica IV/3. 
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is constant in time and it is this value which is to be regarded as the energy 
of the nucleon. In nonrelativistic approximation 


lim (E — mc?) = mv tv : 


C— co 


which is precisely the energy expression of the non-relativistic mechanics. In a 
field depending on time the energy also varies : 


ed (9) 


On taking into account (8) the space components of the equations of motion (7) 
may be written as follows 


d/(E 3 . 
ma | 2 —|1— Seas 
or on taking (9) into account 
2 
d : Al 1 W 
Ip aVv 
mo = = lea + : | (10) 
t eae: t 
mc? 
or in a slightly different form : 
db bp 1 re) | V 
m-—— = —|1 —— rad + — 0 —| mc? ] 1 + — 
dt alle ce Ot oe | or: 


It follows that in the case of small velocities the motion takes place as if the 
particle moved in a field described by a potential 


V* — me? log 


Luare le 


me? 


(This was first pointed out by Plebansky [1].) The equation of motion valid 
for arbitrary energies and velocities is given in (10). 
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3. Motion in a static scalar field 
In the case important in actual practice when in the system of reference 
here employed V is independent of time, the exact equation of motion (10) 
can be greatly simplified. Let us use equation (8) to eliminate the velocity, 


and take into account that owing to (9) E is now constant. Then (10) can be 
reduced to the following simple form : 


|e | [2+ ]eraav, (11) 
a C. 
mo — — grad W (12) 
maior 
me? \2 v2 | 
rot : 


This shows that in a static field the nucleon moves exactly as if the Newton 
force resulting from the potential energy W were acting on it. Starting from 
the Hamilton-Jacobi equation in the approximation E* = E — mc? < mc? this 
potential expression was first derived by Werle [2] in the form: 


y2 
2mc2 


W=V+ (14) 
The present derivation on the one hand supports Werle’s result by starting 
from the equation of motion, on the other abolishes the approximating character 
of that result by generalizing it over arbitrary energies. 

Expression (13) cannot be regarded as a potential in the classical sense 
since it already contains E, one of the constants of motion. Yet, the total energy 
being easily determinable from the experimental conditions in scattering experi- 
ments and by measurements in other cases, the simple equation of motion (12) 
lends itself readily to the discussion of motion, the more as it is identical in 
form with the equation of motion in classical mechanics, and in addition, contains 
no approximation, but describes motion exactly, relativistically. 

From equation of motion (12) one of the integrals of motion is easily 


obtainable in the form 
Ew === > mp2 + WV, (15) 
this being, however, perfectly equivalent to the energy integral expressed in (8). 


3* 
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The calculation is further simplified if it is assumed that y is merely a 
function of distance r from the force centre. Then the equivalent potential W 
also depends on that distance, and the angular momentum is constant. This 
means on the one hand that the motion takes place in the plane, and on the 
other that the absolut value of the angular momentum (which can be 
expressed by polar coordinates) is constant : 


oe == h.= con't. (16) 
dt 
Taking this into account, from (8) 
2 2 2 
dr =el/1- = gee: (17) 
dt E c2 7? 


whence the function r(t) can be readily determined by simple integration. 
Thus the result is that in a scalar static central field the motion of the nucleon 
can be subjected to perfectly exact treatment. 

From the results obtained interesting conclusions can be drawn. Let us 
assume that the potential V is of a purely attractive character, i. e. decreases 
with decreasing of r. Then it is obvious from equation (11) that the accele- 
ration in the region V > — mc? tends towards the centre of attraction. 
However, at the point r, determined by equation 


V (ry) = — me? (18) 


where according to (4) the proper mass vanishes, the acceleration becomes 
zero and within this point for values r < 1%, V < — mc? will then be of oppo- 
site direction. Thus we see that in the case of a potential of perfectly attractive 
character the force centre is likewise surrounded by a repulsive region of a radius 
To, Which is quite independent of the energy or of any other initial condition. 
(A repulsive region exists only if V takes a negative value sufficiently 
great.) 
Let us express the velocity from the energy equation (8) : 


joel) = has: |. (19) 


Hence we can see that at the point r = ry) the velocity of the nucleon reaches 
the velocity of light. This, however, will not cause any difficulty as owing to 
the simultaneous vanishing of the proper mass the value of inertial mass 
M (i — 3/2) remains finite and constant. 
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The nucleon arriving from outside (e. g. in scattering experiments, from 
the infinite), enters the repulsive zone with the velocity of light. In the 
repulsive zone the velocity begins to decrease and if the absolute value of 
the potential becomes great enough it reaches the value : 


V(r) = — m2 — E. (20) 
Then, in the case of radial motion, the velocity of the particle becomes zero, 
i.e, at this point the particle turns back. Accordingly, the motion takes place 
in the following way : the nucleon arriving from the infinite begins to accelerate 


Fig. 1 


rapidly as it approaches the force centre, it crosses the boundary of the repulsive 
sphere with the velocity of light; in the repulsive zone its velocity decreases 
and in radial motion if the potential becomes deep enough, it even reaches 
the value zero. At this point under the influence of acceleration tending out- 
ward the particle turns back and recedes into the infinite. The turning point 
of the nucleon of greater energy lies nearer to the centre of attraction. (See 


fig. 1.) 


4, Examples 


Let us first consider the case when the meson field acting on a moving 
particle is generated by a resting nucleon. Using a scalar meson field we find: 


On the analogy of «classical electron radius» interpreted by the definition 


2 - 
eat a 
Te 
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let us now introduce the concept of the «classical nucleon radius»* r, according 
to the definition 


2 
2 eo = me", (21) 
Ty 
bd pe | 
Considering that g > e, ry is much smaller than the nuclear force range “ , 
wherefore its value is 


o2 2 
s ie ey BS 


As the observed value of g?/ Ac is approximately 1, r,) approximately equals 
the Compton wavelength of the nucleon, i. e. 2. 10~** cm. It will be expedient 
to calculate in 7) units. 

The velocity of the nucleon as the function of r is determined by the formula 
(19). In nucleon-nucleon scattering problems the initial kinetic energy of the 
nucleon is very small as compared to the proper energy, E ~ mc?, thus we can 
write : 


o(r) =<) 1—| a (23) 


mc? r 


(If in the infinite the initial velocity vy = 0, the formula is exact.) The boundary 
of the repulsive zone, where the velocity of light is reached,is at the point 
characterised by the equation 


= Ine" , (24) 


that is, just at the classical nucleon radius, at a distance of 2.1014 cm from the 
nucleon regarded as fixed, thus affording additional justification of the de- 
nomination «classical nucleon radius». (This is about 1/, part of the range of 
the nuclear force.) The radially moving nucleon approaches the fixed nucleon 
over the distance defined by the equation : 


ee 
2 2 J 


= 2mc? 


Og 
| 


Ty 


which, according to 


oA NR u(ro—l 
y= Toe 0 Daz =F, 


* Sat F : 
In fact, to be more consistent, ry should be denominated double nucleon radius. 


; 
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is about one half of the classical nucleon radius. Evidently, in the case of non- 
relativistic velocities two nucleons can never approach each other nearer than 
the half of the classical nucleon radius. 

The motion therefore proceeds as follows. Travelling towards the nucleon 
regarded as fixed the incident nucleon being attracted begins to accelerate at 
about 10712 cm. Within 10-15 cm the rate of acceleration increases inten- 
sely until at about 2.10-14 cm the nucleon reaches the velocity of light. Here 
acceleration must of course stop; on arrival in the repulsive zone the nuc- 
leon loses its velocity and at 10~14 cm turns back and once more traverses the 
region of the velocity of light. Gradually losing its velocity it finally recedes 
into the infinite. We can easily form a notion of motion in time if we consider 
that in the r <1, region e “ <1. Taking (21) into account we can there- 
fore write : 


: 
0) = =e 1 {1 “| : 
dt | re 
Hence we obtain by simple integration : 


1 3 
sea Hote] 
c 


To To 


{t; means the instant of arriving at the turning point r, = r)/2.) From this we 
can see that the time of staying in the repulsive zone is considerably less than 
10-25 sec; during this short time the velocity of the nucleon changes from 
¢ to 0, and again to c. (If the motion is not radial the particle has an initial 
angular momentum and of course reversion sets in sooner, e. g. in p-scattering, 
at about the boundary of the repulsive zone.) 


As a second example let us consider the motion of a nucleon in the field 
of an extensive nucleus. Assuming again a scalar meson field, the potential 
energy in the field of a nucleus regarded as being of constant density, mass 
number A and radius R, will on the basis of (1) be 


3 Ag? 
2 R3 


V (r) = — 


2 (1 + eR) e-#R eal 


ur 
In good approximation it is uR= A" Inthe case of g = 20ethe repulsive zone 


occupies a considerable part of field in the inside of a heavy nucleus. In the 
case of g ~ 5e, however, the repulsive zone does not exist, the equivalent 
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potential (10) does not differ from the potential V(r) and relativistic effects 
do not yield considerable correction. 


* 


In the previous example we assumed the nucleons to be uniformly dis- 
tributed in the nucleus. Let us now examine the motion of a nucleon in a field 


which is generated by other fixed point nucleons. The potential energy. 


of the moving nucleon at the point t in the field of the nucleons fixed at 1; 


(p= 1G 5B) 18 
nh p—# | t—2; | 


it |v — t; | 


According to 3, the nucleon moves as if in a field described in classical mechanics 
by a potential 


= —ty 4 e—2u| 1-1; 
n 2 & e\t t; | § e | AN 


feat R |v —t;| 2me* |t — 1; |? 


gt e# | r—t;] e—#| t—r; | 


= ar (25) 


fFime? |r—r,| |r—y| 


Hence we obtain the important result that the relativistic motion in a static 
meson-field is equivalent to a classical motion in which, in addition to two- 
body forces, there arise three-body-forces: the presence of a third nucleon 
weakens the attraction between the two. (However, this effect is only signi- 
ficant at distances less than 10° 1° cm.) 

Incidentally, a simple calculation can prove that if applied quantum- 
mechanically in the Hartree or Hartree-Fock approximation the relativistic 
effects treated herein lead to three-body forces without any further assumption. 


5. Role of the proper time parameter 


In the basic equation of motion (3) the derivatives of the coordinates 
with respect to the parameter t are involved. The relation of t to the coordinate 
time ¢ is: 


ca ‘pt 
ag os ee (26) 


Hence we see that at the boundary of the repulsive zone at a distance equivalent 
to the classical nucleon radius, where | v| == ao becomes zero, thus at this 
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- point the equation of motion becomes indeterminate. The question arises whether 
this indeterminateness might not be concealing the fact that the nucleon does not 
enter the repulsive zone r < r, at all. To decide it, let us introduce the four- 
momentum p; = Mu;, the components of which may be written as: 


E 
APY Pas Pal, Ys R= —E . 


c 


Although in the case of r = 1, u; becomes divergent, p; does not so because 
in the static field E = const. Let us write the equation of motion (5) in the form 


dp; = F; dt. 


This form does not contain any more singularity, but shows that the nucleon 
traverses the boundary point r =r, with a constant momentum; here dt 
being zero. This means that the nucleon enters the repulsive zone with 
inertial motion without acceleration, as we have assumed on the basis of (11) 
in 3. 

According to (26) the parameter t is to be taken equivalent to the proper 
time of the particle, i. e. the time measured by the clocks in the inertial system 
of reference instantaneously moving with the particle. The differential quo- 
tient of (26) is in general positive, but if the nucleon enters the repulsive zone 
the numerator in the expression of the energy constant in formula (8) previously 
becomes zero. The denominator must likewise become zero. This means that 
(26) is negative in the repulsive zone, t decreases with the increase of t. This 
paradoxical result would appear to show that in the repulsive zone the clock 
of the nucleon goes backwards, which would lead to absurdities in reasoning, 
to the inversion of cause and effect. Let us, however, consider the following: 
we cannot speak of the clock of the accelerating nucleon but only ofthe clock 
of the inertial system of reference, the velocity of which is equal to the 
instantaneous velocity » of the nucleon. According: to the universally valid 
formula of the Lorentz transformation the relation between the time t* 


measured by this clock and between our time t is: 


dt | 


This t* is to be called proper time, which is always the monotonous function 
. . - de* 
of t and generally agrees with the parameter t of the equation of motion ("* /¢, = 
. dr* : . 
= 1); at the boundary of the repulsive zone “* /7, becomes indefinite (%) and 
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within the repulsive zone t* differs from T, (7 la, = — 1). The equation of motion 
can also be written with 1* 

4 (yan) 

dt* dt* Ox; 


, dx i 
This form has the disadvantage that the definition Pj = M as cannot be used 


because when the nucleon passes through the repulsive zone it would result 
in changing the sign of the energy E’ = — pj in static field too, which is always 
i 


constant otherwise. The equation of motion can also be written in the following 
form : 

Ea dp — k OV 

dt* dt* dx; 


pi=|M 
This clearly reveals the cause of the repulsion occuring despite the attractive 


potential: the repulsion is the result of the rest mass becomig negative. 


6. The Hamiltonian 


The Hamiltonian theorem [1] underlying the equation of motion is: 


5| (—me? —V) dr = 0. 


Proceeding to the system time, we obtain the following Lagrangian: 


L =—(mei+V) 


p2 
a. (27) 


Hence the canonical momentum (which is equal to the dynamical momentum): 


Finally setting up the Hamiltonian in the usual way: 


H=py—L= j —— = (me? + V)? + (pe)? (28) 
a ‘ 


c 
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H is equal to the total energy E. The canonical equations naturally lead to 
the equations of motion (10). On the other hand we have seen that with the 
Werle approximation (E* = E — mc? < mc?, but » and V are arbitrary) the 
field equations are equivalent to the equation of motion (12), (14) of non- 
relativistic form, the Hamiltonian of which can be readily obtained : 


Hy, =—-+ BW, Hy =Ey. 
m 


Of course, the two Hamiltonians are not independent of each other, there exists 
the following relation : 
: = met|/ 1+ 2H. 


me? 


thus in the Werle approximation : 
p? 
H=me?+ Hy + maeue Agana pos W|+ Cu 
m | 


Evidently the Hamiltonian leads likewise to the equivalent potential, a re- 
cognition particularly important from the viewpoint of the transition to quantum 
mechanics. (See also appendix.) It is essential to know that this is a concomitant 
to the Werle approximation. If we take the usual non-relativistic approximation 
(b? <c?, V <mce?), then 


r = 
A= me'| (1 eo tet Ss Leet ore 


me? mc 2m 8s 8m c? 2 mick 
? 2 p2V lo" 
L = — mc?}1 — a eo send Py 
mc? Zee 2mc* 8 c? 
(The ... mean terms of order c +.) In this approximation the potential term 


V?/2mc? resulting in repulsion no longer appears. 

(It is interesting to note that the last term is equivalent to one of the terms 
in the Hamiltonian obtainable from the Dirac-equation by means of non- 
relativistic approximation, which has no classical analogue in the case 


of electromagnetic field). 


It is worth while to consider the equations of motion resulting from 
the Lagrangian (27) in non-relativistic approximation. Thereby we may 
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obtain an illustrative picture of the phenomenon of mass defect. Expanding 
the Lagrangian (27) in a power series of 0?/c? and, as usual, keeping only 
the terms c 2, we get 


| 1 b? p?>V 1 b* ; 
V pet A inert 


mc2 2 c2 2mc! 8c? 


te = sd = : s t= —gredV , 
D 
"Qe? 
or 
M’ t = — grad V 
where : 


c? mc2 


which expression gives a demonstrative picture of the mass defects arising 
in the attractive zone of nuclear forces and of the role of kinetic energy. 


7. Motion in a combined electromagnetic and meson field 


Let us consider a particle acted upon by a scalar field and an electro- 
magnetic field (proton). The covariant equation of motion is: 


d Cp e 
= (Mui) = > 2 “4 Fi, Up. 
Tt Oxi c 


Hence we obtain again : 


M =m a as 5 m = const. 
Cc 


We consider the motion in a static field, which is described by a_ time-inde- 


pendent meson-potential @ and an electric potential ©, the vector potential 
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is vanishing. Thus the first three components of the equations of motion are: 


=: BP 

4 ess iy 

¥ oe caer --| 1— 3 &grady —egrad® , (29) 
ey 

whereas the fourth component accounts for the conservation of energy 

4 mc2 dE 

| eee ep er eer (30) 


Hence it is evident that the presence of electromagnetic field does not influence 
the occurence of characteristic relativistic repulsion. The nucleon may accelerate 
until it attains the velocity of light, in the point gy = — mc?, whence it goes 
over to the repulsive zone. In radial motion the turning point is at the distance 
determined by the equation 


gsp—e®=—me’?—E. 


With the aid of the energy formula (30) the equation of motion (29) can also 
be written in the following form : 


ae ie 
S| = =| 1 = 2g grad y — ogra. 
dt 2 c? 


c 


It is rather striking to observe that while the expression of mass proceeding 
that of velocity under the sign of differentiation with respect to time contains 
the mass equivalents of the rest and kinetic energies as well as of the meson 
potential energy, but there does not appear in it the mass equivalent to the 
electric potential energy. In contrast to the other energy forms, in this case 
the Coulomb energy does not lead to mass-defects.* 


* Exact treatment of the problem of mass-defect is possible only within the ranges of 
the relativistic many-body-problem. 
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Pfféets similar to those treated here arise in the case of a nucleon moving 
in a pseudoscalar meson field. We shall revert to them in a subsequent paper, 
where we also wish to treat the problem of radiative reaction. 


re 


We wish to express our thanks to Doz. J. Werle and Doz. J. Plebanski 
(in Warsaw) for valuable discussions initiating the above investigations. We are 
also indebted to Prof. K. F. Novobdtzky and to Prof. L. Jaénossy for valuable help 
in solving important problems of principle which arose during the calculations. 


Appendix 
Equivalent potential in quantum mechanics 
On the basis of the relativistic Hamilton-operator (27) we can readily 


obtain the quantum-mechanical equation of state of a particle of zero spin, 
which moves in the scalar-field : 


o2 
ES ict een (31) 


Avaiee 
5 c2 ot? A2c? 


Let us proceed to the stationary energy-eigenvalue equation. Taking into 
account that in a stationary state y depends on time in the form 


iB} 


a 
p (xyat) = p (xyz) e | (32) 
and introducing the non-relativistic energy according to 
E* = E — mc? (33) 


after simple substitution we get : 


E* ) 2 - 
Be (1+ al v=|— A We 


2mc 


i. 


2mc? | 


Passing over to the Werle-approximation (E* < mc?, » and V are arbitrary) 


- ae 
Ety=|—2— 44m lp, (34) 


‘ 
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_ which is the non-relativistic Schrédingerequation with the Werle potential 
in place of the potential energy. 2 

For the description of a particle (nucleon) of Y4 spin moving in the scalar 
field we can use the Dirac-equation resulting as the «linearisation» of (31) : 


Ox; 


J 0 1 
E +55 (me +V)|p =o. 
, Ac 


The equation can be transformed to non-relativistic form in the well-known 
manner by setting up the following equation : 


Taking into account the commutation relation of the yj-s, and equations 


(32) and (33), we get: 


E* f per isl 


PAT le 


R2 V2 a 
v=| a+\V + = gh Vi ma (34) 


- oa 2mc? 2mc Xj 


Returning to the Werle-approximation this again leads to the classical energy- 
eigenvalue equation of the form of (34) containing the equivalent potential W 
instead of V, apart from a new energy-term depending on the spin. This term, 
the existence of which was first pointed out by Plebanski [1], is in the meson- 
field the analogue of the energy-term leading to the magnetic moment of the 
electron and it results in spin-orbit interaction. 

Thus it is evident that the relativistic term leading to short-range repulsion 
also occurs in quantummechanics. 
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JQBYM)KEHHE HYKJIEOHA B CKAJIAPHOM ME30OHHOM IIOJIE 
r. CAMO u fF. MAPKC 
Pe3ztome 


B | aOoTe paccmoTpeHo KsIaccHY¥ecKOe JBH)KeHHe TOYeYHOrO HYKIeOHa B CKaJLAPHOM 
Me30HHOM losle. YpaBHeHHe ABH)KCHHA B CTATHYECKOM MO0JI€ MODKET ObITh Pelli€CHO COBEPILICHHO 
TOUHO. PacueThI NOATBEPIKAAaOT pesyIbTaTbI Bepse, OTHOCALIMECA K TIPHOIM)KeHHIO MaJIbIX 
gHeprHuit. STH pe3ybTATbI MOKA3bIBaIOT, YTO H B CUyYae YHCTO NPHTATHBaWwUero NOTeHWMana 
(HalIpuMep, B Csly¥ae NoTeHIHana FOKaBa), BOSHUKAIOT OTTAJIKHBAlHe CHAI, CCIM 3HAYeHHE 
NOTCHIMaIbHOM SHEP CTAaHOBUTCA OONbUIe, YeEM 9HEPrHA NOKOA. Takum OOpasom, 6e3 BCA- 
KOrO CHelHasIbHOrO M1, CAMONO.KeEHHA TOMyYaeTCA CYyUIeECTBOBaHHe OTTAJIKMBaWlMxX CHIH, 
paquyc feHcTBUA KOTOPbIX COBMagaeT C KOMMTOHOBCKOH AIMHOH BOJHbI. McceqoBpanne Bpe- 
M€HHOFO XOfla JBHXKCHHA MPHBOAMT K TOMY, YTO YaCTHIa BOKHA NPOXOAUTb Yepe3s rpaHHUy 
OTTAJIKHBaWwWeh OONACTH CO CKOpOCTbHWO cBeTa. B Ciy4ae B3AHMOCHCTBHA MHOTHX HYKJICOHOB 
MOABIAIOTCA IPPEKTLI, CBAZAHHbIC C CHIAMM THIa Tpex Ted. 
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Moffit’s quantum-chemical approximate method was investigated on hand of a simple 
model compound (H,-molecule). The application of Moffitt’s correction to a weaker approxi- 
mate function gives good results, however, the use of a more correct approximate function which 
corresponds to the basic idea of the method, gives a much worse result. We can draw no decisive 
conclusion regarding the method ; the example, however, shows that we must regard the results 
obtained by the method with some caution. On the basis of the calculations it was evident that 
the use of a more correct ionic function gives a much weaker result, thus it is not correct to 
interpret the ground state of the H,-molecule as a superposition of homopolar and ionic states. 


Introduction 


The central problem of quantum chemistry is the search for function 
types from which by use of the variation method we can get functions well 
approximating the wave-function of the molecule. In principle all complete 
orthogonal series of functions which satisfy the prescribed symmetry and boundary 
conditions, are equally suited, because by applying the Ritz-method, and select- 
ing sufficient members, the eigenvalue or eigenfunction can be approximated 
to the desired accuracy. From the practical point of view, however, particularly 
in the case of complicated molecules, the speed of convergence of the series 
is extremely important. In the selection of the suitable form of the function, 
considerations based on physical or chemical experience can give useful hints. 

Recently a new method gained in importance in the theoretical investiga- 
tion of atoms and molecules, in which the wave function is approximated 
by the linear combination of given configurations. This method was applied 
by Jucis [1] and Boys [2] with success to the case of atoms by improving 
the Hartree-Fock method. In the investigation of molecules, a further deve- 
lopment of the molecular orbital method lead to a similar result. [3], [4] In the 
simplest version of the latter the function of the molecule is given as a product 
of one-electron functions (molecular orbital). A further step is achieved if, 
conforming to the wave mechanical Pauli-principle, antisymmetry is required 
in the coordinates of all electrons. The condition can automatically be fulfilled 
if we search for the function of the molecule in form of Slater-determinants ; 
each determinant corresponds to a configuration, that is the electrons are dis- 
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tributed in molecular spin-orbitals. Generally one determinant corresponding 
to a single configuration is not sufficient for a good approximation of the wave- 
function, thus the task is the investigation of the linear combination of different 
configurations. If we search for the molecular orbital in the form of a linear 
combination of atomic orbitals, then the method in this form is totally equi- 
valent to that version of the electron-pair method, in which beside the homo- 
polar function types, functions corresponding to different degree of ionicity 
are also used. 

In the use of the linear-combination of configurations (configuration- 
interaction method) a serious drawback arises from the fact that the number 
of configurations increases extremely rapidly with increasing numbers of electrons 
and we have no criterion as to which configurations are important and which 
can be neglected. The practical use of the method is greatly affected by its weak 
convergence. 

In this situation the initiative of W. Moffitt is very interesting, who 
proposes a new method modifying in some respect the configuration interaction 
method, and improving its convergence. 


Outline of the Moffitt method 


Using thermochemical data Moffitt points out that the total energy of 
the molecule differs but little from the sum of the energies of the atoms con- 
stituting the molecule. Therefore the use of a perturbation calculation is justified, 
which starts from the states of the atoms as being infinitely far from each other 
so that their interaction can be treated as a perturbation. The state of the 
molecule is approximated in the following manner by the linear combination of 
functions constructed from the different states (including the different ionicity) 
of the atoms which are very remote from each other. 

Let us denote the wave function corresponding to the states of the i-th 
atom by @7, pF, .. In the starting state the atoms are infinitely far from each 
other, at the i-th atom there are n; electrons (n; = Z;, because we take into 
account ionic states too, Z; is the charge of the nucleus). Let atom a be in the 
state 4, atom 6 in the state f etc., (a, 8... denote the totality of quantum 


numbers characterizing a state), then the wave-function, representing this 
state of the system is: 


IT -*Scxs D8 (Lat ms) Ong ae Lee t,) shea PRs (1) 


The eigenfunction of the Hamiltonian of the starting state is then: 


He = H,+ Hy +.%.+ He, (2) 
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where H; means the Hamiltonian of the i-th atom. We have 


Be eee We Tee ae (3) 
: H, 0% = w2 Ot (4) 
¥ and 

Wop$ = wa + wh +... + wi. (5) 


‘We call the functions obtained by the antisymmetrization of the product func- 
tions «composite atomic functions» and denote them by yop “71é.In the following 
we denote the single composite atomic functions briefly hyo? 4, Ya, . 

The program outlined above cannot be fulfilled in practice, because, 
apart from the H atom, we do not know the exact atomic wave functions. 
We are forced to use approximate atomic functions, in connection with them, 
however, serious problems arise. The best atomic functions at present available 
(Hartree- Fock method) give very well the values of physical quantities connected 
with the density distribution of the atom, but the term values are reproduced 
less exactly and reliably. The error connected with the latter is 2—5 eV and 
this is of the same order of magnitude as the binding energy of the molecule. 
So if we calculate the elements of the energy matrix (Hp = (Y,HY pdr) by 
means of «composite» functions constructed from the approximate atomic func- 
tions, the same order of deviation is possible in the matrix elements from those 
calculated by means of exact atomic functions, and this fact will strongly 
affect the convergence of the method. 

To overcome this difficulty Moffitt separates the energy matrix into 
atomic and interaction parts. If we use composite functions constructed from 
exact atomic functions, this separation is the following : 


Hag = Has + Vas » (6) 


where H° is the Hamiltonian corresponding to the starting state, V the part 
representing the interaction term V = H — H°. Taking into account that 
H° z= Wp, this can be written as follows : 


Hap = Map We + Vas (7), 
where 


Map = i. WP, aT 


Moffitt remarks that the approximate atomic functions give the density 
distribution fairly well, only the atomic term values cannot be calculated 
correctly. He points out that the correlation effects which are the main cause 
of this deviation are much stronger in the atomic than in the interaction part. 
On this basis he suggests that the Wz values should be given from atomic spectro- 
scopy and the approximate atomic functions used only in the calculation of 
Vag and Mag. (We denote the matrix elements calculated by means of the 


4* 
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approximate atomic functions with ~ to distinguish them from the exact 


values.) 
ive = Map Wes + aes (8) 


The H. a 3 obtained is not Hermitian, thus we must construct the following 


matrix : 
= y geen a : 
Hap eae a (9) 


~The determination of the coefficients of the linear combination can be 
carried out by solving the following secular equations : 


> Ce (Has —E-Maz)=0 . (10) 


The possible energy values are determined by the roots of the secular de- 
derminant : 


| Has —EMap|=0 . (11) 


The separation outlined above is not favourable because it does not permit 
the use of the spatial symmetry of the molecule for the possible factorization 
of the secular equations or the secular determinant. To eliminate this defects 
Moffitt suggests an alternative separation : 


Vas = Has — Map We . (12) 


By Haz is meant the matrix element of the energy matrix calculated by means 
of the approximate functions, Wg is the element of the asymptotic energy 
matrix obtainable from the former if the atoms are infinitely separated. By 
substituting into (8) and (9) and taking into account that Hyg and Maz are 
elements of a Hermitian matrix, we calculate in the secular equations with 
the following energy matrix : 


Has = Hag + [Map (We — Ws) +(W,— Ws)M az). (13) 


This method has the special advantage that it permits an immediate connection 
with the calculations performed up till now, because the My, ... can be made 
to correspond to the starting functions of the method of Heitler-London-Slater- 
Pauling, to that end we must correct the A‘,g matrix relating to them by 
supplementary members. In the latter, besides the elements of the metric matrix 


Bg, 
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(has AB), the Wz foaken from atomic spectroscopy and the We, caleulated by 


means of the approximate functions are necessary. 
Moffitt applied the new method with success to the investigation of the 


_oxygene, ethylene and benzene molecules [7], [8], [9]. In these cases, however, 


it was necessary to use additional simplifications which are independent of the 
basis of the method and which are difficult to control. (Core electrons, the potential 
due to o electrons, etc.) It seems interesting therefore to see the applicability 
of the method in a case, where all the necessary calculations can be performed 
accurately. For this purpose of model-compound the H,-molecule seems to be» 
the best. 


Application of Moffitt’s method to the H,-molecule 


From the point of view of the application of Moffitt’s method the ratio 
of the binding energy to the total energy of the molecule is decisive. In the case 
of the H,-molecule the binding energy is 4,74 eV and the total energy 31,93 eV, 
the ratio is approximately equal to 15%, thus the method could be applied. 
In the case of greater molecules this fraction is smaller, and the use of the per- 
turbation calculation is more justified. This means that the investigation relating 
to the H,-molecule refers to the limit of the applicability of Moffitt’s method. 

According to the basic idea of the method we must construct the starting 
functions from the eigenfunctions of the different states of the atoms thought 
as infinitely separated. The case of the H,-molecule is special insofar as here both 
the ground and the excited states posses exact atomic functions and as regards 
the ionic state we have approximate functions, which give the energy of the 
H -ion very correctly. As the function of the ionic state needs special examina- 
tion, we select first an approximate function for the ionic state which figures 
in the literature and examine the results obtained by the application of Moffitt’s 
correction to this case. 

The simplest possibility is to follow in our investigation the calculation 
of Weinbaum [15], where the state of the molecule is approximated by the 


homopolar + ionic functions : 
y = 61 [tte (1) tp (2) + ta (2) te (1)] + 11 [tea (1) ta (2) + ue (1) uy (2)], (14) 


where 


1 ee 
ua {1)= Be ug ly er es (15) 
|x 
a, 4, mean the distances from nucleus a, b,, 6, the same from nucleus J, relating 
to the first or second electron. We are using in the following atomic units. 


Let us introduce the following notations : 


, ?— nd 
ae 
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P, = ua (1) up (2), Py = uw, (1) ua (2), (16) 
Yq = Ua (1) ua (2), y= (1) up (2). 
y, and are exact eigenfunctions of the Hamiltonian of the starting state 
in which the two atoms are infinitely separated and the 1 electron can be 


found at atom a, the 2 at atom b, or the electrons exchange their places. 
The cor- responding Hamiltonians are as follows : 


Fores ea ree (17) 
2 ay 2 b, 

pea Dea hy: ’ (18) 
2 by 2 Gy 


W, and ~,, however, are not eigenfunctions even in the case of infinite distance 
of the atoms. If for example both electrons are at atom a and atom b infinitely 
remote, then the energy calculated by means of the ~3 = ua(1) ua(2) function 
and with the Hamiltonian 


1 i 1 1 1 
FOO es Ae eee (19) 


2) 2 ay a, Ti. 


deviates by 4,5 eV from the exact energy value of the H ion representing the 
ionic state. 

The application of the twofold version of Moffitt’s method to this case 
means the following : 

a.) Corresponding to each y;the Hamiltonian of the H,-molecule can be 
separated in the following form: H°! + V°'. The explicit form of H®1, H®2 and 
H”° was given further above, H°® can be constructed in an analogous manner. 
We calculate the elements of the energy matrix figuring in the secular equations 
using formulae (8) and (9). W is a diagonal matrix, its elements are as follows : 


W, = W, = 2Ey = — 1 (a. u), W3 = W, = — 0,52641 (a. u).* 


b.) We must use formula (13), the elements of W are given above. The 
elements of the asymptotic energy matrix calculated from the approximate 
atomic functions are the following: W, = W, = —1 (a. a), W3 = W, = —0,375 
(a. u). After determining the elements of the energy matrix, the calculation is 
carried out in the conventional manner. The results of the calculation are 
collected in Fig. 1 and Table 1. In the Figure there are given besides the experi- 


* In the case of the H—ion it is customary to regard as exact value (— 0,52641 a. u) [11] 
that calculated quantum-mechanically by Hylleraas, the experimental value being very uncertain. 


~ 
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mental potential curve the results of the uncorrected calculation of Weinbaum 
and the results obtained by the twofold version of Moffitt’s method. 


| 


Fig. 1. 
TABLE I 
Rm (a. u) D(eV) 
Vern DAUsiit We clara e b.01 6) ere aie 1,677 3,218 
WRHTE Soc5 pe puoues 206 1,400 4,346 
IWIGAINISS | ele oo OF OOUD Lae 1,50 4,940 
1,40 4,740 


liege? pocacooouoe ook aoe 7 


The inspection of the results shows, that the application of Moffitt’s 
method improved the result significantly. Both potential curves approximate 
very well the experimental one. Particularly close is the result obtained with 
version b. The fact that the energy value lies lower than the experimental one, 
does not mean a contradiction to the basic principle of the variation calculus 
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insofar as our calculation is not an exact variation method as we have modified 
single matrix elements. 

The question can be raised, however, whether our good result is a con- 
sequence of the fact that we have taken into account the right energy value 
in the ionic term or whether it is fortuitous. We can decide the question as we 
know an approximate function referring to the H ion, the corresponding energy 
of which deviates only by 0,06 eV from the exact value. This function is the 
following [12]: 


28 —Z(a,;+aQ, 
Ya (1, 2) =— e ees} [1 + yt. +2 (a, — ay)*]. (20) 
: It 


Thus if we use in the ionic part wg (1, 2) for ug(1)u,(2) both conditions to get 
the ionic energy term with good approximation are fulfilled, and we can per- 
form exact variation calculations. 

The calculation can be carried out along conventional lines, the only 
difficulty consisting in the new integrals presenting themselves in great numbers. 
We give in the appendix the formulae for the integrals not figuring hitherto 
in the literature. 

The results of the calculation are collected in Table II. For the purpose 
of comparison we give the results of two other calculations (from these the first 
is due to Scherr [13]) in which the ionic function is substituted by relatively 
simpler ones which give a relatively weaker energy value for the H ion. 


TABLE II 
pa ( 2) Rm (a - u) D(eV) 
(ionic part of the approximate wave function) 
H, molecule 
ce | 
Bg tay +49) 1,491 3,436 I 
M4 
3 
oF (+42) [1 + e749] 1,558 3,254 II 
—- — | 
3 
oF (1*42) [1 + ey 149 + 5 (a, — 04)? i 1,589 3,230 III 
| 


An inspection of the results reveal many interesting features. In spite of the 
use of functions properly approximating the ionic energy, the value obtained 
for the binding energy is but a little better than the results of Weinbaum’s 
calculation, in which he used z = 1 both in the homopolar and in the ionic 
member. 


* As preliminary calculations had shown that the values of the single integrals depend 
very sensitively on the value of the exponent, it was necessary to obtain the value of the expo- 
nents and that of the coefficients c, and c, with higher accuracy than they were known. Table 
III collects these values, 


- er. th ?- , ‘ y r A 
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TABLE III 
i  — 
Ya (1, 2) | = | or Ca : eta! : 

re 
og e*(41 +42) 0,6875 oat — 12,859 
z —z(@\—a9) } 
—é 3) (1c, ria) 0,8257264 0,4933512 — — 13,842 
Gat (a,—a3) | 
ae 1~42) [1 + cy 749+ ¢.(a,—a,)?]| 0,76885 0,318013 0,124718 — 14,291 
13 | 
aie (ME) = ug (1) ua (2) 1 — — —9,52 

' Exact value : — 14,352 


We can compare the results of our calculation with the variation calcula- 
tions using one or two variable parameters, in which the function form is com- 
pletely analogous. These values are shown in Table IV. 


TABLE IV 
z ra 
ug (1) = Vale “y= [ ua (1) up (2) + ug (2) up (1)] + [2a 1) wa (2) + 5 (1) 0 (2) 

z “ R,, (au) D (eV) 
Heitler— London — Sugiura [14] 1 0 1,68 — 3,216 
Wang [15] 1,159 0 1,45 3,75 
Weinbaum [16] £1634 1 | 1,38 3,455 
Weinbaum . [10] 1,193 0,256 1,417 4,026 
Experiment — — | 1,40 4,74 


The variation method using the correct ionic member gives a result a little 
better than the method of Heitler-London-Sugiura, or the simplest version of 
the homopolar+ionic method due to Weinbaum. In both cases the z figuring 
in the exponential is equal to one. 

The result is essentially weaker than in those cases where z (effective nuclear 
charge) figures as a variation parameter. Among the calculations using one 
parameter Wang’s calculation gives a better result. With the given function 
form the best value is obtained by that calculation of Weinbaum, in which z 


and / are variation parameters. 


\ 
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Discussion of the results of calculations; conclusions 


We can draw no unambigous conclusions from the results of the calculations 
as regards the applicability of Moffitt’s method. If we take namely into account 
the original idea of the method, we must use besides the ground and ionic states 
the excited states too. This is especially plausible in the case on the H,-molecule, 
where the energy of that starting state, where the one electron is in the lowest 
atomic orbital and the other is in the first excited orbital, is —0,625 (a.u) 
thus it is lower than the energy of the ionic state —0,52641 (a.u) which we 
have taken into account. We intend to return to the investigation of this problem 
in a subsequent communication. 

Although we can make no decisive statement as regards the method, 
_ still we can draw some conclusions from the calculations made. 

We can state that the results of Moffitt’s method must be regarded with 
some caution. Our result shows that applying the correction to a much worse 
approximate function we get an astonishingly good result, in the case, however, 
where we used a good approximate function conforming to the original idea 
of the method, our result was much weaker. 

The result of this investigation throws some light upon the real role the 
ionic member plays in the approximate function. In the quantum-chemical 
literature the idea appears very often according to which the ground state of 
the H,-molecule can be explained as the superposition of homopolar and 
ionic states. If we compare the results obtained by the three ionic functions, 
we see that the last result was the weaker, the most correct was the H’ function. 
From this we must conclude that the idea mentioned above cannot be correct 
as the function better approximating the ionic state gives the worse result as 
regards the binding energy of the H,-molecule. 

Comparing the data of Table IV the question arises whether it is worth 
while to adopt as starting functions the eigenfunctions of atoms which are 
infinitely separated. The fact that the variation of the effective nuclear charge 
improves the results in both cases so effectively, draws the attention to the 
possibility that in quantum-chemical calculations one must to a greater degree 
than hithertho make use of the effective nuclear charge as variational parameter. 
In the single orbitals figuring in different configurations one has to use different 
effective nuclear charges and although this implies considerable difficulties 
in the evaluation of the integrals, one has to examine, whether this method 
does not lead to a faster convergence in the approximation of the function of 
the molecule. We hope to return to investigations of this kind. 


* 
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discussion of the conclusion drawn from the results of the calculations, to 
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Appendix 


The new integrals occurring in the computations can be divided into three 
groups : 


a) ii [ug (1) ug (2)]? 7”, f(a, ag b, bg) dt, dr, . Notation: "Dx. 
4) § wi (1) uj (1) wi (2) uj (2) r?, f (a, ay by by) dt, dt,. Notation: "EX: 
1c) i) Ua (1) us (2) ug (1) uz (2) r" f (ay ay b, b,) dt, dt,. Notation: "M=: 


In all the three integral types we denote the f (a, a,b, b,) in the following 
manner: the factors figuring in the numerator are written at the right top 


beside the sign of the integral type, and the ones in the denominator are put 
at the bottom. 


(a, = ay)” | 


, ( are denoted briefly by 
(6; — b,)? a 


4 For example : 
¢ b, — A ay 
uz, (1) uj (1) uz (2) uj (2 oyeulbae ue dn Eee 


Tio 


Moreover we use the following notations : 


R Se ht 
a= its. — —(] 2} 9 § = —(]1—:) , 
yarl+z, r=1+R, r=1—R. 


The evaluation of the integrals can be carried out in two ways : 
a.) In the integrals, which contain only three out of the variables a,, 
ay, b,; and by, Ty, can figure as the 4-th integration variable. 


b.) In the other integrals we can use Neumann’s expansion series 


=D ok 2.) Pt (A) Pt (pty) P2 (Up) 005 (21 — $2) 
y=O0 T= 
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where 


Pio 2} Dey =(— 120+ | FEE 


(c+)! 


4, and 4 _ is the bigger or the smaller resp. of the 4, and A,, 4 and m are 
ellipsoidal ‘coordinates. 


a+b a--b 


P’ and Q” mean Legendre polynomials of the first and second kind. Carrying 


out the integration we use the following auxiliary quantities : 


[o} 


% v ak 
W? (m,n, a) eat) e Math 9012). Pia) apae a? —1)2 (= 1) a 
7 tes 


Vv 
Gi (m, a) = [ e "x P* (x) (1 — x2) 2 dx, 


—1 


1 co 
=O ges bc 
ay 3 x" dx, A, (4) = e x" dx, 
=i j 
R 


R 
Cr (aR) = ex" dx , D, (a R) = | ex" dx, 
) 0 


0 


The values of the integrals A,, B,, Gz (m, a) and Wf (n, m, a) can be found 
in tables [17]. 
The new integrals that are necessary for carrying out the calculation 


> 


are as follows: 


owe 
be 1D = raya (He [y* + 2y3 + l6y? + 12y 4+ 112] + 


8y 
+e” [5y + 35]l, y = 2o. 
= 1 @ 0 : 
2. EY = e399 ) 105 Wo (0:0) — 315 Wo (0,2) — 525 Wo (0,4) — 


—175 Wo (1,1) + 1050 Ws (1,3) + 1575 We (2,4) — 1575 W% (3,3) — 
— 189 Wy (0,0) + 630 Wi (0,2) — 525 WY (2,2) + 120 W% (0,0) — 
— 140 W3 (1,1) — 36 W3 (0,0) . 


Ret oa. 7 Te 
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0 


Og a Ome ari 5 Wo (0,1) — 15 Wo (0,3) _—15 Ws (1,2) + 45 Wo (2,3) ++ 
+ 4 W; (0,1) | 
6 
4, 1 E42" Yq = 7 | 10 Wo (0,1) + 30 Wo (1,2) + 9 W2 (0,1) —15 W? (0,3) — 


— 9 Wi (1,2) + 15 W% (2,3) + 4 W2 (0,1) — 3 W? (0,1) — 
eaila) 3.2) . 


re 08 
a 1pAa, as 
Ki 2225200 
— 525 Wo (0,5) + 525 Wo (1,4) + 1050 Wo (2,3) + 1575 We (2,5) — 
— 1575 Wo (3,4) — 378 W4 (0,1) + 630 Wi (0,3) + 630 Wi (1,2) — 


— 1050 W4 (2.3) + 240 W2 (0,1) — 280 W2 (0,3) — 72 W3 (0,1) }. 


} 210 Wo (0,1) — 35 Wo (0,3) — 945 We (1,2) — 


if 8 
6, TER = — 1210 W% (0,1) — 595 W% (0,3) — 525 19 (0,5)— 
z 


— 385 Wo (1,2) + 525 Wo (1,4) + 1050 Wo (2,3) + 1575 Wo (2,5) — 

— 1575 Wo (3,4) — 378 W? (0,1) + 630 Wi (0,3) + 630 Wz (1,2) — 

— 1050 W? (2,3) + 240 W2 (0,1) + 280 W2 0,3) — 560 W2 (1,2) — 
a2 WV (0.1) ; 


3 
7, pan csy, — + ©) _ 499 w3(0,1) + 700 W2 (0,3) -+ 420 W2 (1,2) + 


4.1750 W3 (1,4) — 2100 W% (2,3) + 666 W4 (0,1) — 735 Wi (0,3) — 
— 525 W2 (0,5) — 729 Wi? (1,2) + 315 Wi (1,4) + 840 Wi (2,3) + 
4.525 W? (2,5) — 525 Wi (3,4) — 140 W2 (0,1) + 140 W2 (0,3) + 
+ 240 W2 (1,2) — 280 W3 (1,4) — 36 Ws (0,1) — 36 We (1,2) + 
+ 66 W? (0,1) — 105 W3 (0,3) — 525 W (1,4) + 525 Wy (2,3) — 
— 30 W3 (0,1) + 70 W3 (0,3) +4 W3 (0,1) ( 
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i es 


3, —1Bacosy, — ——<_) _ 490 W% (0,1) — 700 Wo (0,3) + 1820 Wo (1,2) + 


2? 25200 
+. 350 W2 (1,4) — 700 W% (2,3) + 666 W} (0,1) + 105 Wi (0,3) — 
— 525 W? (0,5) — 1569 W% (1,2) + 1155 Wi} (4,1) + 525 Wi (2,5) — 
— 525 W (3,4) — 140 W3 (0,1) + 700 W3 (0,3) — 320 W2 (1,2) + 
+ 280 W2 (1,4) — 560 W2 (2,3) — 36 Ws (0,1) — 36 W3 (1,2) + 
+ 66 W? (0,1) + 315 W; (0,3) — 525 Wy (1,2) — 525 W; (1,4) + 
+ 525 W} (2,3) — 30 W2 (0,1) — 70 Wz (0,3) + 140 We (1,2) + 
+ 4 W; (0,1) | 


1 9 
ie Wa ee eT 1848 W% (0,0) — 9660 W% (0,2) -+ 8085 Wo (0,4) — 


3 176400 

— 3675 Wo (0,6) + 2940 W% (1,1) — 18620 W% (1,3) + 14700 Wo (1,5) — 
—- 16415 We (2,2) — 40425 Wo (2,4) + 11025 Wo (2,6) + 

+ 29400 W9 (3,3) — 44100 Wo (3,5) + 33075 Wo (4,4) — 

— 3780 W} (0,0) + 16884 W? (0,2) — 8820 W7 (0,4) — 5292 Wy (1,1) + 
+ 17640 W} (1,3) — 17640 W? (2,2) + 14700 W; (2,4) — 

— 14700 W} (3,3) + 2860 W2 (0,0) — 9240 W2 (0,2) + 2940 W2 (0,4) + 
+ 3360 We (1,1) — 7840 W (1,3) + 7840 Ws (2,2) — 1120 W3 (0,0) + 
+ 2016 W3 (0,2) — 1008 W3 (1,1) + 192 W% (0,0) | 


1 0° 


Lp ee 1512 Wo (0,0) + 10710 Wo (0,2) + 


2° 529200 | 
+ 5145 Wo (0,4) — 11025 Wo (0,6) + 5586 Wo (1,1) — 39690 W9 (1,3) + 
+ 22050 Wo (1,5) — 6615 Wo (2,2) — 55125 Ws (2,4) + 73500 W4(3,3) — 
+ 33075 Wo (2,6) — 66150 W9 (3,5) + 33075 Ws (4,4) + 2268 W? (0,0) — 
— 17010 W4 (0,2) + 13230 Wy (0,4) — 5670 W? (1,1) + 9450° W2 (1,3)+ 
+ 13230 W4 (1,5) + 29988 W4 (2,2) — 48510 W? (2,4) — 22050 W2(3,5) + 
+ 22050 Wi (4,4) — 780 W2 (0,0) + 2520 W% (0,2) — 2940 W2 (0,4) + 
+ 5160 Ws (1,1) — 10080 W2 (1,3) +° 5880 W2 (3,3) — 168 W2 (0,0) — 
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— 1680 W3 (1,1) + 1512 W3 (2,2) + 192 W$ (0,0) + 384 Weis) 2 
+ 378 W; (0,0) — 4410 W; (0,4) — 882 Wi (1,1) + 8820 W3 (1,3) + 
+ 22050 W; (2,4) — 22050 W? (3,3) — 1890 W (0,2) — 270 W3 (0,0) + 
+ 1260 W2 (0,2) — 1470 W2 (2,2) + 112 Ws (0,0) — 168 W3 (1,1) — 

— 24 W* (0,0) | : 


1 s 
1. *Es, +7E3,= - 


12600 | 440 Wo (0,1) + 70 Ws (0,3) — 525 Wo (0,5) — 
— 2590 Wo (1,2) — 525 Wo (1,4) + 1575 We (1,6) + 3850 W5 (2,3) — 
— 1575 Wo (2,5) — 432 W7 (0,1) + 630 W (0,3) + 2298 WY (1,2) — 
— 420 W{ (1,4) — 2940 W4 (2,3) — 1050 W3 (2,5) + 1050 W (3,4) — 
— 80 Ws (0,1) — 280 W2 (0,3) — 200 Ws (1,2) + 560 W2 (2,3) + 
pate (0.1) 2, (1,2). — 132 WW; (0st) — 210 7, (0,3) 
+ 630 W; (1,2) +1050 W; (1,4) — 1050 W; (2,3) + 60 W, (0,1) — 
— 140 Wz (1,2) — 8 Wz (0,1) . 
2k Ee = = ie | 700 Wo (1,1) — 5600 W% (1,3) + 3150 Wo (1,5) + 
+ 3150 Wo (3,3) — 432 Wi (0,0) + 2100 W7(0,2) — 1050 W4 (0,4) + 
+ 210 W} (2,2) — 2100 W} (2,4) — 760 W2 (1,1) + 560 We (1,3) + 
+ 72 W3 (0,0) — 420 Wj (1,1) + 2100 W; (1,3) + 60 Wz (0,0) — 
— 140 Wz (0,2) . 


The argument of the W (n, m) in 1—12 is always 0 


; co 4 
13. se 
| T=0 =a 
= 41 2(27+1) | 
pis 4 Coe 8 Se (aay 
Bi t* («+ 1)? TD, oul a 


m® (0, a) = G2 (0) [W2 (2,8) + We (1,4) — 2 W? (1,2)] + Gr (2) We (1,2) + 
+ G2 (1) [W2 (2,2) + W (0,4) ~ 2 We (0,2)] + Gr (3), Wr (0.2), 


bg =i LO k . " :- Pi Poae Re! Pe OCR 
; iar * Res k Siow 
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-m® (1, a) = 2 [Gt (1) We (2,3) + Ge (2) We (1,3)], 

m° (2, a) = G2 (0) [2 W2 (0,1) — We (0,3)] — Ge (2) We (0,1) + 
+ G2 (1) [2 W2 (0,0) — Wz (0,2)] — Ge (3) We (0,0), | 

m? (3, a) = — 2 [Gr (1) W2 (1,2) + Ge (2) Wr (1,1) ], 

mt’ (4,4) = — [G2 (0) We (0,1) + Ge (1) Wr (0,0) ], 

mz (0, a) = Gt (0) Wz (1,2) + Gz (1) Wz (0,2), 

mt (2, a) = — [Gz (0) Wz (0,1) + Gz (1) We (0,0) ].- 


The argument of G and W figuring in the m-’s is a. 


162 
14, 1M,, = Ry fe"{a Cy (28) + y?C, (28) — yD, (20) — 4 Do (20) — 


= 15/4y! + et ' AR A, (2a) + y?R3A, (2a) — yR?A, (4a) —4 RA, (4a) | | é 
S2* (0 2p 
Bega | "12 Ga (28) — 2 Cy 26) — 2 Dy (Re) — ry Dy Ce) + 
+ yD, (20) + 2D, (20) —(3 17 + 5) /Ay2)} + e®{27 RAg(2a) + 2R2 A, (2a) — 
— 27 RAy (4a) — 7 yR2A, (4a) — 2R® A, (4a) — yR®A, (4a) | ; 


15.3 MoS 


16. -1M4 = reat [12 C, (28) — 6 C, (2B) + y? C, (28)] — 12 C, (28) + 
+ 6 yC, (28) — y? C,(26) — (5 yr + 9) /4y? — r [3 y D, (20) + 12 Dy (20)] + 
+ 3 y D, (29) + 12 D, (2¢)} 
+ e® \r 12 RA, (2a) — 6 y R? A, (2a) + y? R3 A, (2a)] 
+ 12 R? A, (2a) — 6 y R® A, (2a) + y?R# A, (2a) 
— r[3yR? A, (4a) + 12 RAy (4a)] — [3 y R® A, (4a) + 12 R® A, (4a)] 94 
17, Ma = ne le" {4 r [Ei (28) — Ei ( — 20) + Inz/(1 — 2)] 
— 4 Cy (28) + ry? Cy (28) — y? C, (28) +- 4 Dy (20) — r y Dy (20) 
+ y D, (2) — (15 + 2ry) /4y)} + e® {47 [Ei (— 4a) — Ei ( — 2a)] 
+ 4 RAp (2a) + ry? R? A, (2a) + y? R8 A, (2a) — 4 RA, (4a) 
—ryRA, (4a) — y R24, (4a)* | : 


th | tien ae Ps » ’ a _ ‘ ad 


“Y 
1 
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TT | 
18. 1M = Per Cy (26) + r78G (28) — 4G, 2) 


y 
— 7° C3 (28) — [4 r Do (20) + ry D, (20) — 4D, (2¢) — » D, (2¢)] 
—5(3yr+ 7)/4y?\ + e® {47 RA, (2a) + ry? R? A, (2a) 
+ 4 R? A, (2a) + y?R! A, (2a) — [47 RA, (4a) + 7 y? R? A, (4a) 
+ 4R? A, (4a) + y R® A, (Aa]} 
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UCCIEQOBAHHE S®®EKTHBHOCTH OWHOrO HOBOTO KBAHTOBO-XHMHUE- 
CKOrO MPHBJIWMKEHHOTO METOA 


P. MAYHL 
Pe3rme 


ABTOp HcCeqOBas, KBAHTOBO-xuMuYeCKHH MeTOA B. Modputa WA CayyaA MOMeKYJIbI 
(H,). Mipumenenue nonpasKu Modputa B cayyae He COBCEM TOUHOH MPHONMKeHHOM PyHKWHH 
aeT OTHOCHTeEIbUO XOPOWIM pesybTaT, a IpHMeHeHHe Ooee TOYHOM PyHKUMM, COCTAaBIeH- 
HOM Ha ocHOBe ugen B. Modputa, npuBoguT K yxyfWeHnio pesyibTaTOB. Ha OcHOBe 9THX 
pe3yJIbTaTOB, KOHCYHO, HEIb3A CAeaTb OKOHYATCJIBHbIM BbIBO OTHOCHTEIBHO Mpessaraemoro 
MeTOJa, HO JlaHHbIii MpuMep, BO BCAKOM CJy¥ae, MOKAa3bIBaeT, YTO HeOOXOAMMO NPOABUTD 
NOJDKHY1O OCTOPOXKHOCTh K Pe3yJIbTATaM JaHHOTO MeTOAa. PacueTbI NOKA3AaJIN, YTO TIPHMCHeHHE 
donee TOUHOK HOHHOH PyHKUMH aeT caObIi pesy1bTAT, T. €. HEIb3A PACCMATPHBATb OCHOBHOE 
COCTOAHHE MOJIEKYIIbI BOAOPOMA KaK CyIeplO3HIUMIO TOMCOMONAPHbIX M HOHHBIX COCTOAHHH. 


5 Acta Physica IV/3. 
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In this paper a method for the determination of the sensitivity distribution of photoca- 
thodes is described. The purpose of the investigations was to overcome the difficulties experien- 
ced in reproducing measurements carried out by means of photocells or photomultipliers. The 
variation! of the cathode sensitivity in photomultipliers is determined as a function of location 
-on cathode and of direction of incident light. 


1. Introduction 


In measurements of the intensity of electromagnetic radiation by means 
of photocells or photo-multipliers, the results are often difficult to reproduce. 
This difficulty can be attributed to the following reasons: the temperature- 
dependence of the signal-to-noise ratio of the photocathode, fatigue effects of the 
photocathode [1] and the local dependence of cathode sensitivity. It can also be 
assumed that the photocathode sensitivity is influenced by changes in the direc- 
tion of incident light too. 

In the following we shall at first restrict ourselves to the determination 
of the sensitivity distribution of multiplier photocathodes; further below 
our measurements will also be extended to the dependence of photocathode 
sensitivity distribution on the direction of incidence of light. 

The photocathode sensitivity distribution is determined essentially by 
two components: on the one hand by the local dependence of the sensitivity 
of the cathode surface in relation to the photoelectric effect, (this being due for 
the most part to the inhomogeneous photoelectric properties of the cathode) 
and on the other hand, by the electron optics of the multiplier. In what follows the 
resultant of these two components will be called effective cathode sensitivity. 

The fluctuations occurring in sensitivity distribution play an important 
part especially when the lightcurrent to be measured acts on a comparatively 
small part of the total cathode surface. (The multiplier is used in such a manner 
e. g. in certain astronomical investigations.) In this case the increased difficulty 
arises that the adjustment of the light spot to a definite point of the cathode must 
be reproduced. We note that if the total cathode surface is exposed to incident 
light the role of the sensitivity distribution is less important because the local 


5* 
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differences in sensitivity on the illuminated surface of the photocathode com- 


pensate each other. 


2. Description of the measuring procedure 


For the investigation of sensitivity distribution the following procedure 
has been developed. 

The photocathode is illuminated over a comparatively small area of its 
total surface by a point-like source. The voltage on the anode resistance of the 
multiplier is proportional to the product of the sensitivity shown by the illumin- 
ated cathode point and the incident light intensity. If we keep the intensity 
of the illuminating light source constant the sensitivity on the total 
cathode surface can be determined from point to point. This procedure, 
however, is tedious and its accuracy is not satisfactory. To eliminate this dif- 
ficulty the light spot is moved along a straight line on the cathode surface. 
Connecting the voltage proportional to the displacement of the light spot to 
one of the deflection plate pairs of the cathode-ray tube and the voltage on the 
anode resistance of the multiplier to the other, the photocathode sensitivity 
distribution along a straight line can be directly observed on the screen of the 
cathode-ray tube. During the investigations only constant intensity of the scan- 
ning ray has to be ensured, as we are only interested in the relative sensitiv- 
ity of the individual cathode points. The intensity of the light beam, however, 
may vary within certain limits. 

The block-diagram of the measuring apparatus is shown in Fig. 1. The 
moving point-like source mentioned above has been obtained by exciting the 
screen of cathode-ray tube CRT 1 by an electron beam focussed to about 0,5 
mm diameter and moving the electron beam along the straight line by the 
deflection generator DG connected to one of the deflection plate pairs. 

A moving point-like source developed in the above manner completely 
satisfied our requirements. In addition, its preparation and use is much simpler 
than a corresponding mechanical arrangement would be. 

The image of the light spot moving over the screen of cathode-ray tube 
CRT 1 with an amplitude A is formed on the photocathode of the multiplier 
M with a magnification of N-times. 

For direct plotting of the effective cathode sensitivity curve cathode-ray 
tube CRT 2 is used, one of the deflection plate pairs of which is connected te 
the voltage produced on anode resistor Ry of multiplier M through the wide- 
band amplifier A, while the other plate pair is connected in parallel to the cor- 
responding plate-pair of cathode-ray tube CRT 1. By means of the moving 
point-like light source produced by cathode-ray tube CRT 1 the sensitivity 
of the photocathode is scanned along a straight line. The effective sensitivity 
distribution for the total surface of the photocathode is determined in such 
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a way that by the double potentiometer P a suitable d. c. voltage is connected 
to the other deflection plate pair of cathode-ray tube CRT 1 and on the photo- 
cathode the ray is shifted in a direction perpendicular to the direction of scan- 
ning. Having thus determined the effective sensitivity distribution along several 
parallel straight lines it can also be plotted perpendicular to the direction of 
scanning. This plotting work can be avoided by connecting the deflection 


generator to the other deflection plate pair of tube CRT 1, which enables the 


Fig. 1. CRT1 : 5BP type cathode-ray tube of a persistance time of 10 “sec, with special blue 

screen. CRT 2 ; Double-ray cathode-ray tube (only one of the rays was directly used in the 

measurements). M:931-A type multiplier. L: Mayer Trioplan 1:2,9, f= 7,5 em. Rm = 
= 130 kohms 


sensitivity distribution perpendicular to the former scanning direction to be 
directly observed on oscillograph tube CRT 2. 

During the course of measurements the angle formed by the scanning 
ray and the normal to the surface of the photocathode varies somewhat. Owing 
to this the effective sensitivity distribution may undergo some modification as 
compared to scanning at constant incidence direction. As has been mentioned 
above, our measurements serve, among other purposes, for the determination 
of this effect, a detailed account of which will be given in section 3c. We note 
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here that the variation of the incidence angle does practically not affect the 
results of the measurements. : 

Having assumed the scanning light spot to fall with equal intensity on 
any point of the photocathode it seems to be necessary to investigate the factors 
influencing our assumption in detail. 

First of all, we must point out that the intensity of the light emitted by 
the screen of cathode-ray tube CRT 1 towards the photocathode depends on 
the position of the light spot on the phosphorescent screen. This dependence 
is due partly to the inhomogeneity of the screen and partly to direction- 
dependence of light emission. 

To determine the local dependence of the intensity of the scanning ray 
(intensity-distribution) we have performed a checking measurement described 
briefly in the following. 

The light spot moving on cathode-ray tube CRT 1 with an amplitude A 
is focused on to the photocathode, its image thereby being strongly reduced. 
If the sensitivity of the photocathode can be regarded as constant in the region 
corresponding to the reduced amplitude, the caaulousen on the tube CRT 2 
directly gives the intensity distribution. 

In the arrangement of Fig. 1 with a deflecting voltage corresponding to 
an amplitude A = 100 mm the magnification of the optical system has been 


ih 
selected to be N ~ 40° In this case the amplitude of the displacement of the 


scanning light spot on the photocathode is only 2,5 mm. However, as the sensi- 
tivity of the photocathode cannot be regarded as constant even along a straight 
section 2,5 mm long (see below) we do not employ optical imaging on the cathode, 
decreasing thus the effect of the sensitivity change of the cathode as much as 
possible. The result of the measurement is shown in the photograph taken of 
the screen of oscillograph tube CRT 2 (Fig. 2). In order to determine the relative 
change in intensity distribution the zero level of the intensity is fixed by cover- 
ing up the screen of the oscillograph tube CRT 1 in a width of about 1 mm 
perpendicular to the direction of motion of the light spot. 

As the focusing conditions of the cathode ray tube show a certain depend- 
ence on the anode current the intensity distribution was determined at the 
anode current value used in the determination of the effective cathode sensitivity. 

In measuring the effective cathode sensitivity the value of the anode 
current I, of cathode-ray tube CRT 1 is limited owing to fatigue effects appea- 
ring in consequence of the large currents arising in the last stages of the multi- 
plier [1]. In the measurement of the intensity distribution, owing to the re- 
duction necessary to eliminate the photocathode sensitivity distribution the 
light energy falling on the multiplier is thus for geometrical reasons much 
smaller than in the investigation of the effective cathode sensitivity. Under 
such circumstances stronger amplification of the anode signal of the multiplier 
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is necessary in amplifier A. The latter fact accounts for the higher noise com- 
ponent value in Fig. 2. As can be seen in Fig. 2 the emitted intensity along 
one of the diameters of the screen of tube CRT 1 shows a fluctuation less than 
8%. It must be noted that the latter value sets an upper limit to the fluctuation 
of the intensity distribution because, on the one hand, the photocathode 


[A = 100 mm/ 


Fig. 2. 


Fig. 3. 


sensitivity cannot be regarded as constant even in the 2,5 mm range, and on the 
other, the noise component also contributes to the increase of the measured 
intensity fluctuation. 

The variation of the intensity of the scanning light spot, however, is also 
influenced by the optical system. The effective aperture size of the optical 
system depends on the position P’ of the light spot on see. Thus with 
the angle 6 defined in Fig. 3 the intensity i varies according to 


l = const. cos! 6 . 


In measuring the cathode sensitvity distribution we replaced the optical system 


, 
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employed in the investigation of intensity distribution, the magnification of 

1 ; : : : 
which was N ~ Ae by an optical system of a magnification N ~ ty openly 
tested for the purpose. In the arrangement used for the effective cathode sensi- 
tivity measurement (6)max = 6,5°, the resulting intensity decrease is at most 
2,5%, i.e. it is smaller than the combined effect of the factors investigated in 
the previous control measurement. The size of the fluctuation of the intensity 
due to the above reason and to lens imperfections can be checked with the 
help of the blackening of an emulsion replacing the photocathode. 

We note that the influence of this effect on the results of measurements 
is smaller than the given rate of error as the amplitude A has been determined 
so that even the reduced amplitude VA is larger than the geometrical size of the 
photocathode in the direction of scanning and the maximum intensity decrease 
is obtained on the edges of the photocathode. 

Thus, by means of the above checking measurements we found that the 
assumption for constant intensity underlying the measuring procedure is satis- 
fied with an inaccuracy of about 8%. 

In what follows we shall briefly summarize the aspects decisive for the 
selection of the data of the electronical arrangement shown in Fig. 1. 

The anode circuit time constant of the multiplier and the repetition 
frequency f of the deflection generator were determined on the basis of the 
accuracy requirements of the measurement and so were the frequency and phase 
characteristics of the employed amplifier. 

Let us assume that the effective cathode sensitivity rapidly descreases 
in the direction of scanning over a section of length s, correspondingly the voltage 
on the anode resistance of the multiplier decreases exponentially with a time 
constant t = Ry. Cm, where Cy is the total capacity of the anode point to 
earth. Knowing the constant deflecting velocity v of the point-like source on 
the cathode-ray screen the scanning time corresponding to the distance s is 


To attain the purpose of our measurements we require 
{SG 5 


In the choice of fit must be kept in mind that owing to the focusing con- 
ditions of the cathode-ray tube the mean diameter of the scanning ray on the 
photocathode is about 0,18 mm, i. e. it does not seem suitable to select for s 
a value smaller than 0,2 mm, for the given magnification. Accordingly, we have 


chosen f = 215: Hz, Ry = 130 kohms at a capacity of about Cy < 40 pF. 


” 
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In order to check the complete electrical system in the arrangement 
shown in Fig. 1 |and at a magnification of N ~ 6 the scanning light spot on 


the screen of cathode-ray tube CRT 1 is periodically extinguished by means 
of an impulse generator connected to the grid of the cathode-ray tube and 
synchronized by the deflection generator DG. The image then appearing on 
the screen of cathode-ray tube CRT 2 is shown in Fig. 4. The latter method 


Fig. 4. 


also serves to check to what extent our conditions for the neglecting of 
the persistence of the screen in connection with the point-like light source are 
fulfilled. 

As the oscillogram shows the electrical circuit satisfies our requirements, 
as an intensity decrease corresponding to s = 0,6 mm is reproduced almost 
without distortion on the screen of cathode-ray tube CRT 2. 

The procedure described above essentially satisfies the Pe Tae emicnis Oy 
a measuring. procedure suitable for the determination of cathode sensitivity 
distribution. These requirements are 

a.) short measuring time . 

b.) recording of the result of measurements by pceimirabion; a 

Our procedure had to be developed in view of the time EOC DIAS: of 
commercially available multipliers and light sources, EUS RE Me if, for a) 
reason, a multiplier of definite cathode sensitivity distribution is needed, owing 
to differences in manufacture many multipliers would have to be tested for 


selection. 
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3. Results of measurements 


The investigation carried out by the measuring procedure described in 
section 2 can be divided into three groups. 

a.) Scanning of the photocathode of the 931-A type multiplier in parallel 
with the axis of the multiplier (see fig. 5a, direction a). 


a 
A 
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Fg focusing grid 
1 dynode 17. 


nt light | 
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Fig. 5a and 5b 
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b.) Scanning of the photocathode of the 931-A type multiplier perpendicu- 
lar to the axis of the multiplier (see fig. 5b, direction 6). 


c.) Investigation of the dependence of the effective cathode sensitivity 
of the multiplier on the direction of the incident light. 


Fig. 6a 


Fig. 6b 


We note that the investigation in c.) is needed not only in connection 
with the use of the multiplier but it is also of importance in checking the accuracy 
of the a.) and b.) type measurements. . 

To a.) The measurements have been carried out im an arrangement cor- 
responding to Fig. 1 ata magnification of N = 0,31 by connerune the Volizze 
supplied by generator DG to a deflection parallel to the axis of the multiplier 
on cathode-ray tube CRT1. In this case the image 1s formed on photocathode P,, 
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of the multiplier (see fig. 5b) in such a way that the light spot is imaged on a 
clouded glass replacing the photocathode. The effective cathode sensitivity 
curves for two different 931-A type multipliers are shown in the photograph 
taken of oscillograph tube CRT2 (fig. 6a and 6b resp.), where the scatter of the 
sensitivity distributions of the two multipliers can be well distinguished. 
From several measurements we have obtained for the size h of the sensitive 
surface of the photocathode parallel to the axis of the multiplier on the average 
20 mm, while the geometrical size is h) ~ 25 mm. 


Fig. 7. 


As can be seen in the oscillograms the effective cathode sensitivity as 
a function of longitudinal size shows strong fluctuation. The effect of the eloceras 
static focusing grid placed in front of the cathode of the multiplier (see fig. 5a) 
can be well distinguished in the figures, the number of minima of the cathode 
sensitivity curve being in good agreement with the number of grid windings. 

Tob.) The measurements have been carried out with an arrangement sinnilae 
to that used in a.) with the difference that the roles of the vertical and horizontal 
deflection plates of cathode-ray tube CRT1 are interchanged and for the ampli- 
tude the value 45 mm is chosen. We have obtained for the projection of the sensitive 
cathode surface size perpendicular to the symmetry axis of the multiplier t = 
9 mm (see fig. 5h) whereas the corresponding geometrical size is tf) ~ 9 mm. 
The measurements resulted in the oscillograms of oscillograph tube CRT2 
presented in Fig. 7 at various a values (see fig. 52). The variation of the curve 
shapes with a comparatively small change in the value a is due to the effect 


of the grid windings forming a small angle with the scanning direction. In the 
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latter measurement the focusing is aimed at approximately the bisecting 
point of the cathode surface size perpendicular to the symmetry axis of the 
multiplier. The geometrical size t, of the cathode projection being comparatively 
small, greater defocusing does not occur. We note that for the estimation of 
the error thus committed the b.) type sensitivity curve can be plotted from 
the a.) type measurements. 

To c.) The measurement is carried out with an arrangement similar to 
that in b.) with the difference that the multiplier is turned about its symmetry 


Fig. 8. 


axis by an angle a (see fig. 5B). The photographs taken of the effective cathode 
sensitivity curve at various values of a are shown in Fig. 8. With the aid of 
these photographs the angle-dependence of the effective sensitivity for a definite 
point of the cathode surface can be approximately determined (Fig. 9). The 
determination of the accurate angle-dependence of the effective sensitivity is 
made very difficult by the electrostatic focusing grid. 

As in measurement b.) the scanning ray and the given direction of incident 
light form a maximum angle of 4° the variation of the effective aperture size 
influences our measurements by 1%. 

In measurements a.) the direction of the scanning ray and the normal 
to the surface of the cathode form a maximum angle of 6,5°. In the checking 
of our measurements of the intensity distribution the effect of this varying 
angle of incident light was also investigated. It can be seen that the accuracy 
of our measurements was not influenced by the change in the direction of 


incident light. 


4 
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By means of the measuring method described above, in addition to the 
cathode sensitivity investigation of the various types of photomultipliers and 
photocells, other optical adjustment and control problems can as well be solved: 

My thanks are due to Prof. L. Jdnossy for valuable direction during the 
course of these investigations as well as to the Industrial Research Institute 
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for Telecommunication and especially to Academician G. Szigeti for having 
prepared the short-persistence screen cathode-ray tube used in our measurements. 
I am also indebted to J. Kmetyé for help in building the measuring equipment 
and to B. Lendvay and Maria Tassaly for excellent measuring work. 
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ONPEHEJIEHHE PACIPEQEJIEHHA UYBCTBHUTEJIBHOCTH KATOWHOM 
MNOBEPXHOCTH ®OTOYMHO}KUNTEIIA 


*K. HAPAM 
Pe3tome 


ABTOPOM OBL. paspaboTaH MeTOA, NOAXOAAMMH WA ONpexeneHuA pacnpenenenusa 
YYBCTBUTCIbHOCTH (OTOKATOAHBIX NOBepxHocTeH. Llemb1o paspaboTKH 9TOrO MeTOZa ObINO 
BbIACHCHHE TPyHOCTeH, BOSHHKAIOWIUX B CBASH C NOBTOPACMOCTh M3MepeHHit (OTOSTeMeHTOM 
MIM POTOYMHOXKUTeIeM. C NOMOLMIbIO 9TOTO MeTOAa ABTOP ONpeseH H3MeHeEHHe UYBCTBUTECJIb- 
HOCTH B 3ABHCHMOCTH OT MeCTa HM OT HalipaBIleHHA Maarojero cBeta. 
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DIE STATISTISCHE THEORIE 


DES ATOMKERNS 
III. Teil 


Von 
P. GOMBAS, E. MAGORI, B. MOLNAR und FE. SZABO 


PHYSIKALISCHES INSTITUT DER UNIVERSITAT FUR TECHNISCHE WISSENSCHAFTEN, BUDAPEST 


(Eingegangen: 1. X. 1954) 


Wahrend in den vorangehenden Arbeiten die statistische Theorie des Atomkerns fiir die 


_ skalare Yukawasche Nucleonen-Wechselwirkungsenergie vom Typ J = — ye—|t—t'\/ro/) r—r’| 


entwickelt wurde, wird hier die statistische Theorie des Atomkerns fir die skalaren Nucleonen- 
Wechselwjrkungsenergien J = — ee—| t—t'|/"0, sowie J = —ee—(|t—t'|/")* weiter ausgebaut. 
Es zeigt sich, dass die mit dem empirischen Befund am besten iibereinstimmenden Resultate mit 
der Wechselwirkung J = — ee—|t—t'|/ro erzielt werden. 


In zwei vorangehenden Arbeiten wurde von einem der Verfasser unter 
Voraussetzung einer skalaren Yukawaschen Wechselwirkungsenergie zwischen 
den Nucleonen von der Form 


pee (1) 


das statistische Modell des Atomkerns entwickelt! und die Energie des Kerns, 
sowie die Nucleonendichte im Kern mit Hilfe des Ritzschen Variationsverfahrens 
in erster, zweiter und dritter Naherung bestimmt. Diese Naherungen unter- 
scheiden sich darin, dass in der ersten Niaherung eine Gausssche Dichtever- 
teilung mit einem einzelnen Variationsparameter, in zweiter Naherung ein 
verallgemeinerter Dichteverlauf mit zwei Variationsparametern angesetzt 
wurde, der sich vom Gaussschen Verlauf beliebig unterscheiden kann, und in 
dritter Naiherung die Neutronen- und Protonendichte voneinander unabhangig 
variiert wurden. Im Falle der genannten Yukawaschen Wechselwirkung fiihrt 
die zweite und dritte Naherung praktisch zu keinerlei Anderung gegeniiber 
der ersten und fir die Nucleonenverteilung ergibt sich nicht nur fir die leichten 
sondern auch fiir die schwersten Kerne eine Gausssche Verteilung, die sich nur 
darin unterscheiden, dass die Dichteverteilung der schweren Kerne mit wach- 
sendem r bedeutend langsamer abklingt als die der leichten. 

In der vorliegenden Arbeit soll nun mit demselben Verfahren wie in den 
genannten Arbeiten die Energie und Nucleonendichte fiir die skalaren Wechsel- 


1 P, Gombés, Acta Phys. Hung. 1, 329, 1952 im folgenden als I zitiert und Acta Phys. 
Hung. 2, 223, 1952 im folgenden als II zitiert. Man vgl. weiterhin P. Gombds, Nature (London) 
170, 362, 1952; 171, 979, 1953; Ann. d. Phys. (6) 10, 253, 1952 ; 12, 156, 1953. — Alle Bezeich- 
nungen der vorliegenden Arbeit sind dieselben wie in I und II. 
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wirkungsenergien von der Form 
Viera coma © (2) 
J = — ae te We) (3) 


berechnet werden, wo « hier eine Konstante von der Dimension einer Energie 
bezeichnet; rp) ist wieder eine Konstante von der Dimension einer Lange, tiber 
die wir hier noch verfiigen kénnten, wir wahlen sie jedoch wieder so wie in I, 
d.h. wir setzen sie mit der durch 2 dividierten Comptonschen Wellenlange 
der z-Mesonen gleich. In der vorliegenden Arbeit sollen nun die Berechnungen 
fiir diese Wechselwirkungsenergien in erster und zweiter Naherung auf genau 
dieselbe Weise wie in den genannten vorangehenden Arbeiten durchgefihrt 
werden. Alle Bezeichnungen sind dieselben wie dort, so dass wir uns durchweg 
auf diese beziehen und uns sehr kurz fassen kénnen. 

Der einzige Unterschied gegeniiber I besteht darin, dass jetzt in der 
Dichte der Austauschenergie statt der Funktion f(@,, w,) von I die folgenden 
Funktionen stehen. 


Fur den Fall der Wechselwirkung (2) 


“0 | don op — [1 +3 (0% + 08) If In [1 + (On + O)4 — 


4x3 } 


—In [1 + (@, — wp)? | } + 4 (03 + w}) aretg (Wn + @),) — 


f (@n,@p) rae 9 


-—4 (w3 — 3) erctg (MW, — Wp) | 
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und fiir die Wechselwirkung (3) 


Ona Me 


f (@n,@p) = = | Ps (w3 | w3) @D 
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bezeichnet ; w, und wp, haben dieselbe Bedeutung wie in I. 


DVR ei Di a 


DIE STATISTISCHE THEORIE DES ATOMKERNS III. TEIL So 269 


\ 


Die Berechnung der Kernenergien wurde in zweiter Naherung fir die 
Kerne A — 16, 80 und 200 durchgefiihrt. Die Resultate sind zusammen mit 
dem Wert von / in den Tabellen 1 und 2 dargestellt. Fiir die Kernladungszahl 
Z wurde hierbei ganz analog zu I jeweils der Wert gesetzt, den wir in erster 
Naherung aus der Minimumsforderung der Energie erhielten. Die Resultate 
sind zum Vergleich mit den in I erhaltenen und den experimentellen in Fig. 1 
auch graphisch dargestellt. 

Wie zu sehen ist, erhalt man fiir die Kernenergien fiir die beiden Wechsel- 
wirkungen (2) und (3) mit dem experimentellen Befund etwas besser iber- 


Fig. 1. Die auf ein Teilchen entfallende Kernenergie E,,/A in zweiter Naiherung (in MeV-Ein 
heiten) als Funktion von A, fiir die Kerne, die fiir ein vorgegebenes A die tiefste Energie besitzen 
O mit der Wechselwirkung (1) berechnet, 
[] mit der Wechselwirkung (2) berechnet, 
VY mit der Wechselwirkung (3) berechnet, 
@ empirische Werte. 
Die mit der Wechselwirkung (3) berechneten Energien (VV ) fiir 4 = 80 und 200 fallen praktisch 
mit den Stellen [] zusammen und wurden der besseren Ubersicht halber nicht eingezeichnet. 


einstimmende Werte als mit der Yukawaschen Wechselwirkung. Besonders 
hervorzuheben ist, dass fiir die Wechselwirkung (2) die Energie der leichten 
Kerne bedeutend tiefer liegt als im Falle der Yukawaschen Wechselwirkung, 
wodurch sich fiir diese Kerne insbesondere im Falle der Wechselwirkung (2) 
eine bedeutgnd bessere Ubereinstimmung mit den experimentellen Werten 
ergibt als fiir die Yukawasche Wechselwirkung. 

Die Lage der Isobaren mit tiefster Energie im (A,Z)-Diagramm 
erweist sich fir die Wechselwirkung (3) und insbesondere fiir die 
Wechselwirkung (2) ebenfalls als etwas besser als mit der Yukawaschen . 


Wechselwirkung (1). 


6 Acta Physica IV/3. 
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TABELLE 1 


° . ° ° \ . PP fiir einig e 
n und Werte der Variationsparameter in erster und zweiter Naherung 
aoe Kerne im Falle der Wechselwirkungsenergie (2). (A= 5,36). 
Alle Energien in MeV-Einheiten. 


A 16 80 200 
Z 8 36 82 
Erste Niherung (y = 0) co 3,78 3,90 3,85 es 
4,73 5,85 5,95 
Zweite Naherung { He 0,54 14 1.8 
E.JA erste Naherung 7,857 8,177 7,064 
ool A ¢ sweite Naherung 7,908 8.651 | 7,735 
TABELLE 2 


Kernenergien und Werte der Variationsparameter in erster und zweiter Naherung fiir einige 
Kerne im Falle der Wechselwirkungsenergie (3). (A = 6,13). 
Alle Energien in MeV-Einheiten. 


A 16 80 200 
Z 8 36 83 
Erste Niherung (y= 0) cy 5 5.2 SIS 
Zweite Naherung } - ie tae = 
Rod erste Naherung 6,987 7,955 6,660 
wl A ) zweite Naherung 7,007 8,728 7,807 


Gemiss dem Ansatz fiir die Nucleonendichte [man vgl. I (135)] ergibt 
sich fir die Nucleonendichte fiir die beiden Wechselwirkungsenergien (2) und (3) 
in erster Naherung eine Gausssche Verteilung und zwar durchweg fiir alle Kerne. 
Im Gegensatz zu I und II erfahrt jedoch die Nucleonendichte der mittelschweren 
und schweren Kerne in zweiter Naherung eine wesentliche Anderung und zeigt 
fiir diese Kerne mehr keinen Gaussschen Verlauf, wahrend die Nucleonendichte 
der leichten Kerne auch fiir diese Wechselwirkungen einen Gaussschen, Verlauf 
aufweist. Fiir mittelschwere Kerne ergibt sich im Kerninneren ein konstanter 
Dichteverlauf, wihrend der Dichteverlauf der schweren Kerne im Kernmittel- 
punkt eine Mulde aufweist (man vgl. Fig. 2), die darauf zuriickzufihren ist, 
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dass die Coulombsche Protonenabstossung fiir die schweren Kerne so betracht- 
lich wird, dass die Protonen durch diese Abstossung an den Kernrand gedringt 
werden und zufolge der starken Anziehungskrafte zwischen Protonen und 
Neutronen auch einen Bruchteil der Neutronen mit sich ziehen. Als Folge 
dessen entsteht am Kernrand eine geringe Erhéhung und im Reramictelpankt 
eine geringe Verminderung der Nucleonendichte. 


ae oes eS 6 


r— 


Fig. 2. Verlauf der Nucleonendichte 2 = On-+ @p in zweiter Naherung fiir die Isobaren tiefster 
Energie mit der Massenzahl 4A = 16,80 und 200 im Falle der Wechselwirkung (2). 


im r, Einheiten. 


Mit den Dichteverteilungen der zweiten Niherung wurde auf dieselbe 
Weise wie in I und II mit dem Zusammenhang (210) in I (7 = 4) der Kern- 
radius R berechnet. Die Resultate befinden sich fiir die zweite Naherung in den 
Tabellen 3 und 4. 


TABELLE 3 


Kernradius in der zweiten Naherung fiir einige Kerne im Falle der Wechselwirkungsenergie (2). 
R in 7,-Einheiten. 


Theoretisch R 
A a 
vA R in zweiter Naherung mayiiecs 
16 8 3,9 2,65 
80 36 6 4,52 
200 82 8 6,14 
TABELLE 4 


Kernradien in der zweiten Naherung fiir einige Kerne im Falle der Wechselwirkungsenergie (3), 
R in r, Einheiten. 


el 


Theoretisch R 
A 5 ; _ empirisch 
Z R in zweiter Naherung 
16 8 2,8 2,65 
80 36 4,5 4,52 
200 83 6,5 6,14 


Der Unterschied zwischen dem Verlauf der Nucleonendichten, die sich 
fiir die Yukawasche Wechselwirkung einerseits und fir die Wechselwirkungen 
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(2) und (3) anderseits ergeben, ist darauf zurickzufihren, dass die Yukawasche 
Wechselwirkung fiir verschwindende gegenseitige Entfernung zweier Nucleonen 
singular wird, wahrend dies bei den Wechselwirkungsenergien (2) und (3) 
nicht der Fall ist. Zufolge dieses Unterschiedes ergibt sich auch zwischen den 
in all’ diesen Arbeiten zugrunde gelegten Austauschwechselwirkungsenergien der 
Nucleonen ein betrachtlicher Unterschied, der fiir den Fall einer konstanten 
Nucleonendichte und N —Z = A/2 aus Fig. 1 in I ersichtlich ist. Aus dieser 
Figur ist klar zu sehen, dass fiir diesen Fall die auf Grund der Wechselwir- 
kungen (2) und (3) berechneten Austauschwechselwirkungsenergien der Nucleo- 
nen fiir grosse Nucleonendichten in eine Konstante iibergehen, d.h. dass 
diese Wechselwirkungen einen Absattigungscharakter besitzen. Im Gegensatz 
hierzu verlauft die auf Grund der Yukawaschen Wechselwirkung der Nucleonen 
berechnete Austauschwechselwirkungsenergie fir grosse Nucleonendichten 
wesentlich anders, und zwar geht sie in eine steil abfallende Gerade tber, 
woraus ersichtlich ist, dass fiir die Yukawasche Wechselwirkung keine vollstan- 


dige Absattigung der Krafte besteht. 


CTATHCTHYUECKAAH TEOPHA ATOMHOLO APA III. 


ll. FOMBAILL, 9. MATOPH, B. MOJIBHAP, E. CABO 
Pe3tome 


B npeabiyyumx padorax Oba pasBHTa CTaTHCTHYeCKaA TeOPHA aTOMHOTO AApa ANA 
CKaJIAPHO'O MOTeHUMasIa HYKICOHHOTO B3anMOelCTBHA THNa KOKaBa J = — ye—|t—v'|/7o/jy—r’\. 
B sToH Ke padoTe Oba MOCTpoOeHa CTATHCTHYeCKaA TeEOPpHA AIA CKaNApHOrO NOTeHWMana 
B3auMomeicTBuA THNa J = — ée—|t—v'l/to, a TakKoKe AIA TuMa J = — ee—(\t—vt’|/to)*, OKasbi- 
BaeTCA, YTO Pe3yMbTATHI, HaMOOMee COBMAaaWUIHe C SKCMePHMeHTOM, MOMy4awTCA B Ciyyae 
creqyroujero MOTeHWMaa B3saimMOsetcTBuA | = — ee—|t—r'l/To, 


UBER DEN 4/-ZUSTAND DES O;-MOLEKULS 


c Von 
A. BUDO und I. KOVACS 


(Eingegangen: 8. XI. 1954) 


SPEKTROSKOPISCHE ABTEILUNG DES PHYSIKALISCHEN ZENTRALFORSCHUNGSINSTITUTES 
DER UNG. AKADEMIE DER WISSENSCHAFTEN, BUDAPEST 


Die am 4//,-Zustand des O, Molekiils beobachtete Multiplettaufspaltung, die als Funk- 
tion der. Rotationsquantenzahl von den auf Grund der van Vleckschen Theorie abgeleiteten 
Formeln abweicht, kann durch eine Stérung des 4J/-Zustandes von zwei 2J/-Zustanden in hin- 
reichend! guter Ubereinstimmung mit den Messergebnissen gedeutet werden. Aus der Theorie 
dieser St6rung lassen sich die angenaherte Lage und die Multiplettkonstante des einen, bisher 
noch nicht beobachteten */7-Termes angeben. 


Unter den Elektronenzustanden der zweiatomigen Molekile kennt man 
bisher nur einen +/I-Zustand, dessen Feinstruktur genauer untersucht wurde; 
das ist der 4//,-Zustand des OF Molekiils, das untere Elektronenniveau der 
von Nevin analysierten ersten negativen (42, — 4I/,,) Banden [1]. Die experi- 
mentell gefundene Multiplettaufspaltung dieses *//-Termes zeigt in ihrer Ab- 
hangigkeit von der Rotationsquantenzahl J keine gute Ubereinstimmung mit 
den Formeln, die Brandt von der van Vleckschen Theorie ausgehend abgeleitet 
hat [2], [3], sondern es ergibt sich eine Abweichung von einigen cm 1, im 
wesentlichen von solcher Art, dass die mittleren zwei Komponenten der eine 
von den dusseren naher kommen, als es theoretisch zu erwarten ware. 

Diese Abweichungen sind nicht dadurch bedingt, dass die Brandtschen 
Formeln die ersten Glieder einer Reihenentwicklung sind ; rechnet man namlich 
auf Grund der erweiterten Formeln 
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"ls 
W,= Fad) = Bl +) — 2 | + AI +H) 5 — 28] + 
ppt), ¢ b= MEE), 
nt ST+ Int 4tU+0% 
3 23 26 |7/: 
W.= FU) =B|JU+N— S| t IAD + Oty 
BP See 1 ieee 
ATT 1) ee a0 
: A 6Y (Y +4) 
= gh ge = (¥ ley ie 
[wely, = (Y — 4), yo (Yi -7, - oy Fares 


ist und B die Rotationskonstante, A die Multiplettkonstante bedeuten], oder 
lést man die entsprechende Sakulargleichung fiir einzelne J-Werte numerisch 
auf, so bedeutet dies nur ganz unwesentliche Anderungen. Die Diskrepanz 
zwischen den experimentellen und theoretischen Werten ist umso bemerkens- 
werter, weil sich die van Vlecksche Theorie im Falle der //-Terme niedrigerer 
Multiplizitat recht gut bewahrt hat. 

Fir die Erklarung der Abweichung kann die Wechselwirkung der 4//-Term- 
komponenten mit anderen Molekiilzustanden in Betracht kommen. Durch die 
Beriicksichtigung der Wechselwirkung mit einem +2’-Zustand gelang es den 
Verfassern, die A-Aufspaltungen beidem in Rede stehenden 4//-Term in Uber- 
einstimmung mit den Messungen zu deuten [4], und dieselbe Rechnungen ergeben 
auch Verschiebungen in der gegenseitigen Lage der Termkomponenten, doch 
sind diese Verschiebungen gegen die genannten Abweichungen viel zu klein 
(sie bleiben unter 0,1 cm™~1). Es zeigte sich, dass die Wechselwirkungen mit 
einem anderen 4//- und mit einem ‘4-Zustand ebenfalls zu vernachlassighar 
kleinen Verschiebungen fiihren. 

Es gibt aber Wechselwirkungen auch zwischen Termen verschiedener 
Multiplizitat, die durch den Operator 2a; (I; S;) der Bahn-Spin-Kopplung der 
einzelnen Elektronen vermittelt werden [5], so dass zur Deutung des fraglichen 
Problems die Wechselwirkungen des 4//-Zustandes mit 2-, 2/7- und 24-Termen 
in Frage kommen. Unter Beachtung der entsprechenden Auswahlregeln [s. die 
spatere Gl. (24)] ergibt sich, dass nur der *JJ-Term eine Verschiebung der 
beschriebenen Art in der gegenseitigen Lage der 4/7 -Komponenten bewirken 
kann. Das wird auch dadurch unterstiitzt, dass im Termschema des O3 auch 


mehrere *//-Terme zu finden sind, die verhaltnismiassig nahe an dem betrachteten 
4/T-Zustand liegen. 


F 
"4 
d 
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Zur Erlauterung und mit Riicksicht auf das Folgende sind die Term- _ 
schemata von O, und OF — nach Mulliken [6], aber durch die neueren Ergebnissen 
erginzt — in Abb. 1 dargestellt. 

Da nach dem Obigen die beobachtete anomale Verschiebung der 4/J-Term- 
komponenten durch ihre Wechselwirkung mit einem 2/7-Term (bzw. mit mehreren 
*J1-Termen) verursacht werden kann, wollen wir im folgenden diese Wechsel- 
wirkung niher untersuchen. 

Die Lésung der Wellengleichung des zweiatomigen Molekiils hat bekannt- 
lich unter Vernachlassigung gewisser (die Separation hindernden) Glieder in der 
Gleichung die Form 


YW — DRu; (2) 


dabei ist ® die Lésung der Wellengleichung 
[H — V(r)] d=0 (3) 


des Zweizentrensystems, u und R sind die Eigenfunktionen der Wellengleichung 
des symmetrischen Kreisels bzw. des anharmonischen, Oszillators. Wir wollen 
den Teil des Hamilton-Operators H in (3), in dem Spinglieder nicht vorkommen, 
mit H,, die entsprechenden Higenfunktionen mit ®, bezeichnen : 


[Hy — V (r)] %=9, (4) 
und den erwahnten Operator 
H, —— Pet) (I; S;) (5) 


der Bahn-Spin Wechselwirkung als Stérung betrachten; die im allgemeinen 
viel kleineren gemischten Glieder mit |; $; und die Spin-Spin-Wechselwirkung 
kénnen hier ausser acht bleiben. 

Die erste Aufgabe ist die Bestimmung der Eigenfunktionen ®) derjenigen 
Elektronenzustande von O3, die bei unserem Problem eine Rolle spielen. Das 
sind finf //-Zustande, die wir im folgenden der Kiirze halber mit 


i lies Aly ATT (6) 


bezeichnen. Den Grundzustand *//* — nach der Mullikenschen Bezeichnung 
...(wx)* (vc) X24, — kann man nach Mulliken aus dem Grundzustand 
. . . (wx) 4(vm)? 8X7 des O, entstanden denken, indem man ein Elektron aus der 
Ausseren v-Schale entfernt. Ahnlich gehen die iibrigen vier Zusténde von (6) 
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Elektronenvolt 


1,0 1,5 2,0 25 
r(cm x 1078) 


Abb. 1. Ein Teil von den Termschemata (Potentialkurven) des O, und O%-Molekiils. Die 
Zustinde X?II,, A®JIp und der gestrichelt eingezeichnete, bisher noch nicht beobachtete 2/7 
Term von O, entsprechen der Reihe nach den spater im Text mit 2/77, 277°, 277° hezeichneten 

Zustinden, 


ee ee 
gS 
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aus der (wz) *(vz)?-Konfiguration des O, durch Entfernen eines [je mach den 
einzelnen Zustanden verschiedenen | Elektrons aus der w-Schale hervor — alle 
vier Zustinde haben also die Konfiguration . . . (w2)3(vm)? —, und zwar sollen 
die bekannten 4J/- und 2/7°-Terme in diesem Sinne dem Term 3d", von Q,, die 
unbekannten *//°- und *//*-Terme den metastabilen 12} baw. 1A, von O, ent- 
sprechen (vgl. auch Abb. 1). 

Bei der Bestimmung von ®, gehen wir von den Eigenfunktionen der 
einzelnen Elektronen aus, die in Zylinderkoordinaten 0, z, g die Form 


p (0) e'“?a, baw. p (z, 0) e'“?B haben. a und f sind die zu den Werten + - 
bzw. — & der Spinquantenzahl o gehérenden Spineigenfunktionen, fiir 7-Elektro- 


nen ist A = + 1; die Funktion p beziehe sich auf die w-Schale, die entsprechende 
Funktion fiir die v-Schale sei r. Die vier méglichen Eigenfunktionen der Elektro- 
nen der w-Schale sind also pe'?a, pe~‘?a, pe'?f, pe ‘8, diejenigen der v-Schale: 
re'?a usw. 

Sieht man zunachst in (4) von der elektrostatischen Wechselwirkung 
zwischen den Elektronen ab, so lassen sich bekanntlich Eigenfunktionen fir 
die einzelnen Konfigurationen in Determinantenform angeben. Betrachten wir 


1 
7b.die Zustande mit A= 24,—= + 1,2 = 20; = soe (also mit 2 = A + Y= 
3 
—— a so ist die entsprechende Eigenfunktion fir die Konfiguration (wz) 4(v7) 


(unter Fortlassen der Elektronen der inneren abgeschlossenen Schalen) : 


pie*'a, pre Pa, pze'a, pye't*a, pe "a, 


[Pie By Peep, pre "B, pre ™B, pse Bs 
1 


1f*e+-- j 
Dy = Pi ar Phe, pre PG, Ps¢ "Gs Ps ea, Ds en a, ~ Ys! A l(a, 20,84) le a 


V5 


| Pi eB, Pe ane P3 oma Pe mB, Ps one , (7) 


slag ips ips 
lr,ePa, r,e'Pa,. r,e'P'a, rea, re Pas 


[In der leicht verstandlichen abgekiirzten Bezeichnung [ | der Determinante 
bezieht sich also die erste runde Klammer auf die w-, die zweite auf die v-Schale. | 


Fir die Konfiguration [ebenfalls mit 4=+1,2=-4- | ergeben sich vier 


Eigenfunktionen, und zwar mit den abgekiirzten Bezeichnungen die folgenden : 


1 nS ee 
yal Alle|| ’ 


2 [teidlec, m= Plot. Alead. 


Na Bs l(c, 6 Bs az) (a, é,)I. D, = 


ie (8) 


~ Ys 
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Die finf zueinander orthogonalen Eigenfunktionen D; sind, wie man leicht nach- 
weist, zugleich Eigenfunktionen des Operators L; = 2’ I,; mit dem Eigenwert 


1 
A = 1 und die des Operators S, = Ys,; mit dem Eigenwert Y = cs Die Eigen- 


funktionen der Zusténde (6) miissen aber ausserdem auch Eigenfunktionen 
des Operators S$? = Ss? +S} + S? mit den entsprechenden, die Multiplizitat 
der Terme bestimmenden, Eigenwerten S(S + 1) sein. Die Eigenfunktionen, 
die auch der letzteren Forderung geniigen, die also schon den Zustanden (6) 
zugeordnet werden, kénnen, seien mit B; bezeichnet. Nun ergibt die Anwendung 
des Operators S? auf die D; folgendes : 


S*Dy= i Dy = [5 + 1}Do, also (HIf,) By = Do: 


|: 


S*D, =A [7 D, + Ds + Dy 


7 
SD, =i (D, + {Dat Da) 


3),=i [D, + D, +73,| : 
S?D,=??- “Dy also” (@U73,): By =D. 


Wahrend also D) und D, die Eigenfunktionen der Zustande *//;,, und 27s, dar- 

stellen (wie erwahnt, ohne Beriicksichtigung der Wechselwirkung der Elektronen), 

hat man, die Eigenfunktionen von 4//, 277° und 2J7° als Linearkombinationen 
3 


> ¢D; aus der Forderung 
1 


3 3 
S? Se; Di = 779A S cj Di [mit 4= S(S+ 1)] (10) 
1 1 


zu bestimmen. Dies fiihrt zu den Gleichungen 


c 
pjatata=o, 


(11) 
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Die Wurzeln der Sakulardeterminante und die entsprechenden S-Werte sind : 


ee SO 8 sere gay 
| 4 2 4 2 


Die Sreel A, au der die Lésung cf!) = cf) = cel gehort, entspricht also 
dem */7,,-Zustand, mithin ist seine normierte Eigenfunktion : 


Cia 


a 


(D, + D, + Ds). (12) 


Fir die’ zur Doppelwurzel 2? = 4 gehérenden Zustande 2/7” und 2/7° folgt nur: 


ec — = (i + c?), cS) = — (c+ eS), d. h. [mit ce? ) = ¢;, cf or ci]: 
(?/7;,,) : By = ¢, D, + ce Dy— (ce, + Cy) Ds ; (13) 
(7U2;),): Bs = c, Dy + ¢2D,— (c; + ¢2) Dy. (14) 


Die Eigenfunktionen dieser Zustande bleiben also nach diesem Verfahren 
noch unbestimmt, und die Beriicksichtigung der elektrostatischen Wechsel- 
wirkung der Elektronen, durch die die Entartung aufgehoben wire, lasst sich 
quantitativ kaum durchfiithren. Wenn wir aber die Mullikensche Zuordnung 
zu Hilfe nehmen [nach der der 2//’-Zustand in dem erwahnten Sinne dem 3Z- 
Zustand von O,, und 2J/° dem 1X-Zustand des O, entspricht] und beachten, 

Halse ee 
dass zur (vz)?-Schale beim 12’ (O,) nur die Bestandteile acd (8 a) gehéren, 
die nach (8) nur in D, und D, auftreten, so kénnen wir annehmen, dass die 
Higenfunktion des dem 12 (O,) entsprechenden ?//;,-Zustandes nur D,- und 
D,- Anteile und keinen D,- Anteil haben kann. Setzt man dementsprechend 
in (14) c; + cp = 0, d.h. c, = — cj, so ist die normierte Eigenfunktion 


Die Eigenfunktion B, von 7//° ergibt sich daraus, dass B, und B, zueinander 
orthogonal sein miissen ; das fiihrt zu c; = c, und nach Normierung zu 


1 
Det Dy 2 D),)'- 


6 


( a1T3,) 3 B, sce 


f 
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Ahnlich kann man auch die B-Funktionen der anderen Komponenten 
der in (6) stehenden Terme erhalten. Wir fassen die Resultate zusammen : 


(:II2,): By =Dy: (212,): Bo = Do: 
(417;,,) ‘ B, = D,; 5 (427_1),) 5 B; — dD. 3 

7 1 oa i 4 
(*22s),) > B= ae + D, a Ds) > (*10,,,) : B,= ya + D2 + D3) 3 


(15) 
vA 1 vd - 4 
Clay Bs == oe (D,+ D,—2D,); (7/12): Bz =—(Di+ Dz2—2D3) ; 


/6 /6 


Vs! 


und die Funktionen D; gehen aus den entsprechenden D; durch «Spiegelung» 
der Spinfunktionen — d. h. Ersetzen jedes a, durch $; und jedes 6, durch a, — 
hervor. Samtliche B-Funktionen sind orthogonal zueinander und normiert. 

Die Funktionen B; wiirden, wie erwahnt, im Falle der Vernachlassigung 
der gegenseitigen Wechselwirkung der Elektronen die Eigenfunktionen ®, 
von (4) darstellen. Wir betrachten die B; als Naiherungseigenfunktionen und 
berechnen mit ihnen die Matrixelemente des Operators H, = 2a; (I; S;), und 
zwar zuerst die Diagonalelemente fiir die Termkomponenten 


(@Hi,): By = . (D, — D,) : (:U1,): Bs “5 Ce ay. 
(II4,): B, = Dy: CIt,): Bi = Di; 
dabei ist 
Doone (16) 


‘ / o77o d 
*TTay5 thle, 5 Aida, Pion lage 


Die Diagonalelemente haben bekanntlich die Form AA, wo die Multiplett- 
konstanten A bei mehreren Termen aus Messungen bekannt sind. Auf Grund 
der folgenden Rechnungen werden eben diese Konstanten mit gewissen, Mittel- 
werten der im H p stehenden a; in Beziehung gesetzt, wodurch auch die Bestim- 
mung der nichtdiagonalen Matrixelemente von H, miglich wird. Fiir die Diago- 


nalelemente ergeben sich auf Grund der bekannten Eigenschaften der Operatoren 
I; und s; die folgenden Resultate. 


a ee 
“ ad 
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1) ay] (A= sls 1, += a Tye A?ya AX’ == cee a i Bj a: we Bo dis 


Nei 
ae ~;, lle Ey eal 2a; 1; $2 ibe aylle af eee eee P; % pi — 


* im * eek * 4 
— Pr % Pk + Pj 2 Pi + Pm om Pm + Tp An tn) = ets aa: +ap+ap -+-a,) = 


i] (17) 


a : 
— G,; 


die Bezeichnung ¥'P( ) bedeutet dabei das Summieren des in der ( ) stehenden 


-Ausdruckes iiber alle Permutationen von (iklmn) = (12345), ferner ist 


@, == 7, AnTn = Fe? i) T* (On, Zn) @ (On; Zn) ¥ (Ons Zn) On dOn den (18) 


und entsprechend a,. Unter Verwendung des experimentellen Wertes von 


A2zya [8] besagt also (17) : 

Aza G7 == 195 em}. (19) 
Zum gleichen Resultat gelangt man natiirlich durch die Betrachtung der UT, 
Komponente mit Hilfe der Eigenfunktion B,’ in (15). 


2) 4//-Zustand. Hier (und bei den folgenden zwei Zustanden) lasst sich 
die Berechnung der Diagonalelemente nach (15) auf die Integrale 


‘ 1 
| Dj Hp D, dt =F (4 + 2a), {|D3H,D, dt peed a 2a;\; 
(20) 
{ D3 Hp D, dt = > Cai | D7H p Dy dt = 0 (i,k =1,2,3;i £ k) 


zurickfiihren. Die erhaltene Beziehung nebst dem beobachteten Wert von 
Az [8] lautet : 
1 
17 oe in ap = — 48 cm}, also a, = + 144 cm +. (21) 
3) Bei dem 2//’-Zustand ergibt sich 


1 
Asgr => 4 (22) 


d.h. nach (21) sollte Az,» = +48 cm ! sein, wahrend der experimentelle 
Wert 8,2 cra—! betragt [8]. 
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4) und 5) Bei den unbekannten 7/7" und 2/7"-Zustanden fiihren die Rech- 


nungen zu 
Aarne Se See Apo A2yd —— —_ ap ° (23) 


Unter den nichtdiagonalen Matrixelementen von H, = 2a; (I;S;) kénnen 
bekanntlich nur diejenigen von Null verschieden sein, die den Ubergiangen 


AQ = 0, 4A 0, 41, 48 = 0,451 (24) 


entsprechen. Die Rechnung ergibt, dass in unserem Fall auch einige von diesen 
Elementen verschwinden ; die von Null verschiedenen sind : 


2 
Hy (*Hs),, *11%,,) = — Hp (*Hy,, 21.) = — s ap = — 67,88, 


UY 
Hp (4115, *1T5/,) Sas Hp (*1L,),, *1T5),) = 2 Se aes a 159,22 ? (25) 


- Hy (UT3,,, 21T5),) = — Hp Ui, *2Tt,,) = ice 1 112.588 


V3 


Die Zahlenwerte der Matrixelemente (25), die fiir das Spatere von Wichtig- 
keit sind, kénnen aber nicht hinreichend zuverlassig erscheinen, da die Ab- 
weichung des theoretischen Wertes (48 cm +) von Aa, in (22) vom beobachteten 
Wert (8,2 cm~!) darauf hinweist, dass die ohne Beriicksichtigung der gegen- 
seitigen elektrostatischen Wechselwirkung der Elektronen erhaltenen Eigen- 
funktionen (15) modifiziert werden miissen. Die elektrostatische Wechsel- 
wirkung bewirkt, dass die richtigen Eigenfunktionen (bei einigen von den 
betrachteten Zustanden) die Linearkombinationen der Funktionen B; werden. 
Da bekanntlich die elektrostatische Wechselwirkung schon in viel einfacheren 
Fallen nicht genau in Rechnung zu tragen ist, suchen wir die Koeffizienten 
der Linearkombinationen aus der Forderung zu gewinnen, dass der theoretische 
Wert von Az,» mit dem beobachteten iibereinstimme. 

Die elektrostatische Wechselwirkung ergibt bekanntlich Matrixelemente 
nur zwischen solchen Zustaénden bzw. Funktionen B; einer gegebenen Elektronen- 
konfiguration, zu denen dieselben Quantenzahlen A, Y und S gehéren, in unserem 
Fall also zwischen, den Termkomponenten 2/7 ie 277. wes 277, (und zwischen he 
21T;,, J1%,). Daher hat die Eigenfunktion ®, von Ty, die Form [vgl. (15)] =: 
®° (2IT:),) = S2By + SspBs + SapBa, und ahnliche Form haben auch ©, (2/I<,) 
und @, (217%, ). (Soo, S3o,.-. sind die Koeffizienten der Linearkombination.). 
Zur Bestimmung von Sop, S35,... reicht offenbar die Kenntnis der erwahnten 
Konstante Az,» allein nicht aus, und es liegt keine andere verwendbare Grisse- 
vor. Wahlt man nun nur einen von den beiden Zustanden TT), und tae so: 
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_ fihrt — wie leicht zu zeigen ist — die Wechselwirkung des 2/7’ aly ‘Zustandes 
mit dem 2/7 vl zu dem der Beobachtung widersprechenden Wert Az,5> ae 


Mit Ricksicht hierauf nehmen wir eine Wechselwirkung nur zwischen den 
277° und 2/7°-Zustinden an, d.h. wir setzen: 


Dy (217) = So, B, + S3B,, 
(26) 
DATION Sa Bye Sa By 


Die Koeffizienten S lassen sich aus den puthogonalitatss und Normierungs- 


bedingungen (S2pSec + SspSgc = 0, ne + Sn el Sz: + See = 1) und aus der 
erwahnten Bedingung 


e : 1 Bere 
{Do (2/7:,,) H,® (?/2;),) dt = Seo re «| —- Sso rE 3 ap Se 282S:5| 


= As AS = 3 As p= 4,1 (27) 
bestimmen ; wir geben gleich ihre Zahlenwerte an: 
Sap = Sao —-0,99633, S3p = — Soe = — 0,08557 . | (28) 
Mit diesen S-Werten hat man ausser (26) auch die Eigenfunktionen 
®, (21T;,.) = Sop Bz + SspBs 3 By (210i,,) = Soc Bz + Sac Bs (29) 
zu bilden, die ®,-Eigenfunktionen der iibrigen Zustinde sind nach unserer 
Annahme mit den entsprechenden B-Funktionen identisch. 


Mit den Eigenfunktionen ®, gerechnet, bleiben die Werte von A»,», A 
Aza ungeandert, fiir Az,» und A». gilt (27) bzw. 


ieee 


Aye = 8, ts “>| Bese eas eA Ste Sec a = 1) (30) 


und fiir die nichtdiagonalen Matrixelemente ergibt sich statt (25) : 


H, (4Us,, 2173+) = — Hy (Ty, IT) = Soy |? or —S»!? a, =o = — 81,26, 
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H, (*Hs,,, *1T5),) = — Hp (41Ty, Ty) = Sse |/2 a, — Soc 12 ap =o=+ 152,82 , 


i 
Hy (2112),, 2118,,) = — Hp (217%,, 2174,) = Sas Se bes «| — S3, S3-ap + (31) 


te 


y3 


+ (S25 Sse + $3; S2e) = a, =t = + 127,30, 


d.h. die kleinen B,- bzw. B,-Anteile der Eigenfunktionen (26) rufen nur ver- 
haltnismassig kleine Anderungen in den Werten der nichtdiagonalen Matrix- 
elemente vor. 

Wir gehen nun von dem bisher betrachteten Zweizentrensystem zu dem 
schwingenden und rotierenden Molekiil itiber. Als seine ungestérten Eigen- 
funktionen wahlen wir mit Riicksicht auf (2)—(4) die Funktionen 


Y,—@,Ru, (32) 


und betrachten H, + H, als Stérung, wo der Operator H, die bei der Separation 
der vollstandigen Wellengleichung vernachlassigten Glieder [2] reprasentiert. 
Dann ergibt sich, dass die hier untersuchte Wechselwirkung der ‘//-Termkompo- 
nenten mit den *//-Termen in erster Naherung durch die Sakulargleichung 


aT ol Wie ee 0 0. ru 

417, Neo Hae a Weg Tae 0 cape Hise 

‘Tl, | O He. Ws—W Hi, 0 Hive 20 H3. 

AT ikO 0. «HE 70 O O Onan 

Raleianpa wid pls, She wom a FSS ale ht Sane Shane ete aig pelt wien = ee en = —— O (33) 

=H Hees BL Hs, O OCW Hea se 

217? | O O Hes O | Hes We-W: 0 He, 

aT 0 Hea 0. “He, Oo Sy he 

a IO O He, 402 ee) Hee: Hay eee 


beschrieben wird ; hier bedeuten : 


19 


3 
W? = Bay e— 2) + 5 Aer WY (IT; I) 4 Bao eae Arb ; 


th b < “4, 
y ‘ ; : 
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5 | 
W3 = By fet | tg fis B= 9 OT 0) + Bane 2] — Ay ; 
: 4 


13 1 Ww? 
W3 = Bey et | an P= (UEs TF) + But [2 — a Ray: 
2 


5 3 
Wo = By at a ou Weel; 277°) + Beye Ft yl> Arye; 
(34) 
: sali well 
' Ht= B|/ 3|2—12 3 35 = Bry» | x2: 
ers: 
3, = 2B |x — = i= Bie |) vag} 
| iS 
Hi = Bry | 3 z+}: Hg, = — Heg=7; 
H}s = — H35 =e ; Ay, =—Hxy=c; HY = Hp, 


und x= J (J+ 1). Vonden Abstanden » (4//, 2/7 4 und » (4//, JI) der ent- 
sprechenden Terme ist der erstere —6900 cm™!; der letztere sowie Bzyc sind 
nicht bekannt. Die fiir die entsprechenden Terme mit A = — 1 giiltige Deter- 
minantengleichung ist mit (33) identisch. 

Die Diagonalelemente von (33) sind die Energiewerte, die sich aus der 
separierten Wellengleichung des Molekiils ergeben und die dem Hundschen Fall 
a) entsprechen ; bei ihnen ist die infolge der Rotation auftretende gegenseitige 
Wechselwirkung der Multiplettkomponenten nicht bericksichtigt. Diesen 
Wechselwirkungen entsprechen die Matrixelemente innerhalb der durch ge- 
strichelte Linien getrennten Quadrate ; die iibrigen Matrixelemente ergeben 
sich dadurch, dass der Operator H, =)" a; (I; $,) in die bei der Separation der 
Wellengleichung vernachlassigten Glieder aufgenommen wurde, bzw. dass wir 
von den Eigenfunktionen (32) ausgegangen sind. 

Transformiert man die Matrix H° von (33) mit einer solchen unitaéren 
Matrix S, dass die 4/7 und 2/7-Anteile von H° Diagonalform annehmen, so stellen 


in der transformierten Matrix 
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W, O O Hoot Wilts Hy, Ay, 

O W, O O Ay, gg Hy, Fg, 
Ovsrexc 0W 2 oat Hoo. Hea Hew 
Oe Oe eee Wee Hes ee tae Be Soe (35) 
Hy Hig Ha Hie Ware 0 eee 

He = Hes. “Hest Hage Oe Wiens | Tee Bee 

ig hey Az, A, HH, i ee W, O 


Hy, gy Figg Fig, Ag; Hyg 0 W, 


die Diagonalelemente W,, W,, W;, W, gerade die Brandtschen Energieaus- 
driicke (1), die Diagonalelemente Ws, W, und W.,, W, die bekannten Hill-van 
Vleckschen Formeln [2] fiir die. *[T’- und “II°-Terme dar, ferner sind 


Hus = Sax Ss5 H8; + Ssx Ses W963 Her = Sex S77 H3; + Sax Ser H3,; 


Hye = Sox Ss¢ H8; + Sax See H953 Hus = Sox Hrs HE; + Sax Sse H3; (36) 
(i= 1 Se4) 


Hiz = Ss; S77 S87 + Se: Sar HS; His = Ss; S78 HS; +) Sor Sas |Hes 
(l= 5,6). 
Fir die Elemente von § gilt [4], [9] 


Sik = [sf + sox + sx see sax)” Sik (i, ke Lhe 2, 3, 4) ° 
3 5\1*| 0 1}? 
xe -[5/7- 5] +2] \or3—wa) ors —Wi)—3[7 +] |; 
2 
sox = —(W9 — Wi) |ors—Ws) (73 — Ww) 3/545 (37) 


vx = 2|(J -+] [s+ alk (Wm —w)W—-Wy, 


ae 37S |[e+s]] [y+ 5 |ore—wi) ; 


Von (S55 S56, Sg5> Seg) und (S77, Sze, Sgz, Sgg) geniigt es fur das Folgende nur 
soviel zu wissen, dass sie je eine unitare Teilmatrix bilden ; die ubrigen Elemente 


von § sind Null. 
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; 


3 


a 


Wendet man nun auf die Matrix (35) die iblichen Formeln der Stérungs- 
theorie von nichentarteten Systemen an und bedenkt, dass die gegen die Multi- 
plettaufspaltung grossen Termdifferenzen » (4/7, 2/7”) und v (4/7, ?/I°) fir alle 
Termkomponenten als gleich genommen werden kénnen, so erhalt man fir 
die Energiewerte der ‘//-Komponenten : 


W, = W,+0(Si+S3,) (k= 1,2, 3,4) (38) 


_-mait 


o2 s oe 
-» (417,217) sy (417, 217°) | 


(39) 


Nach dem Obigen ruft also die Wechselwirkung des 4//-Zustandes mit den 
*/I-Termenin dem gegenseitigen Abstand der *//-Komponenten die Abweichungen 


AW; — AWis = (W,, — W;) — (Wi — Wi) = a (S3;, + S3, — S3; — S3;) (40) 


von den auf Grund der Brandtschen Formeln (1) erwarteten Werten vor. Die 
beobachteten Abweichungen sind fiir den Schwingungszustand v = 0 des 
477-Termes, nach der Tabelle X X der Nevinschen Arbeit, in Abb. 2 durch Kreise 
dargestellt ; die ausgezogenen Kurven zeigen die theoretischen Werte, die 
aus (40) mit der Konstante a = — 3,712 cm [nach numerischer Berechnung 
von S;,] gewonnen wurden. Man sieht, dass die beobachteten Abweichungen 
durch die untersuchte Wechselwirkung im wesentlichen erklart werden kénnen. 

Wir bemerken, dass die Giltigkeit der Resultate von (38) bzw. (40) von 
den unter (14) und bei (26) gemachten Annahmen nicht berihrt wird; ohne 
diese wiirden nur in der Konstante a in (39) auch andere Gréssen auftreten. 
Den mit diesen Annahmen erhaltenen Ausdruck von a kénnen wir aber dazu 
beniitzen, um gewisse Aussagen iiber die Lage des bisher nicht beobachteten 
2]7’-Termes zu machen. 

Zunichst folgt aus der Richtung der beobachteten Abweichungen [die 
mittleren Komponenten von 4// kommen der tiefsten naher], dass a<0 ist. Da 
nach der Analyse der 2J7’-?/7*-Banden [7] der *JI ’-Term iiber dem 4// liegt [d.h. 
» (477, 211°) < Oist], gilt nach (39): » (4/7, 2JI°) <0. Bei der Bestimmung von 
y (4IT, 21I°) beachten wir, dass das v=0-Niveau des 477-Termes von allen Schwin- 
gungsniveaus von 2J’und 2/I° gestért werden kann. Dementsprechend schreibt 


man genauer statt des ersten Gliedes von a: 


SY [JR50(r) Re dr fuTudtoe]? — SV PLS RORe dr]? _ P 


(41) 
v’=0 y (411,03 27°,) v=0 w(All =o; 2]7°,) 


9 


7* 


288 A. BUDO und I. KOVACS 


AW,2 -AW,2 


4W,3-4W23 


= men gto bre. 2° 191 °b_6° 2-2 0101S8-1 926% 3-19 — 


OF Se 0 \ 15) $20 eS rao eas SIes 

i i n die in den gegenseitigenAbstanden der 4J7-Komponenten beobach- 
ullgh eee nie nach der Peres Theorie erwarteten Werten, abs der Reihe nach 
die aus der Tabelle XX der Nevinschen Arbeit [1] entnommenen Differenzen AF,;,0ps (J)— 
—AF 4s, calc (J), 4F 39; obs (J) —A F g9,calc (J), 4F5,,0bs (J)—A F,,,calc (J). Die ausgezogenen Kurven 
stellen die durch die Stérung des*//-Zustandes von zwei?J/-Termen verursachten, aus der theore- 
tischen Formel (40) berechneten Verschiebungen 4 W’,,(J)—4 W,.(J), 4 W3'(J)—A W23(J) und 

AW,4(J)—AW4(J) dar. 


wo Ry und R,, die Schwingungseigenfunktionen von 4/7 (v” = 0) baw. 2/I’(v’) 
bedeuten und 


Le eet ta eee oe ree 
y (*H 2032115) Seve (tll. Fe We E +5 aa 


het 
+ We Xe 


1 \2 
x sae o) 42 
w+ (42) 


v, (412; 2T’) = — 6895 cm™}, w, = 1035,69 em}, wx, = 10,39 em}, w, = 
= 900 cm, w,x, = 13,4 em71 


sind [8]. 

“4 sich die Kernabstinde der ‘JI und 2J1°-Zustiande in der Gleichgewichts- 
lage (r, = 1,3813 bzw. 1,409 A) sowie ihre Schwingungskonstanten w, nur 
relativ wenig voneinander unterscheiden, kann man annehmen, dass das Schwin- 
gungsniveau v’’ = 0 des ‘II hauptsachlich vom Niveau v’= Odes 2/7” gestért 
wird. Im Falle von kleinen Schwingungsquantenzahlen v’ kénnen die Eigen- 
funktionen Ry, R, des anharmonischen Oszillators angenihert durch diejenigen 
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des harmonischen Oszillators ersetzt werden. Die mit den letzteren nach den 
Hutchissonschen Formeln [10]* berechneten Integrale {Ro R,-dr fahren zu 


St ye ROR dre 
v=o 9 (417,03 T°.) 


0,0001 + 0,0003 + 0,0003) = — 0,96. (43) 


= pw 


— (0,8598 + 0,0942 + 0,0010 + 


Was nun den unbekannten 2/7 ’-Term betrifft, kénnen wir annehmen, dass bei 
ihm der Kernabstand und die Schwingungskonstante nicht wesentlich von 
denjenigen der 4/7- und 2/7°-Terme abweichen, da ja alle drei Terme zu derselben 
‘Elektronenkonfiguration gehéren. Das bedeutet, dass das Niveau v’’ = 0 des 
*1T-Termes im wesentlichen ebenfalls nur vom Niveau v’ = 0 des 2/I° gestért 
wird. Der Abstand des letzteren von dem ‘4//,_9 Niveau ergibt sich aus (41) 
und nach dem Obigen zu 


2 
Ati pone) Ia oe ig. ~~ — 8500cem—?. (44) 


a— 


Nach unseren Uberlegungen sollte also der bisher nicht beobachtete 7/7° Zustand 
(in der Abb. 1. gestrichelt eingezeichnet) um etwa 8500 cm! iiber dem 4/I,_ 9- 
Niveau, oder [nach (42) ] um etwa 1700 cm ! iiber dem bekannten ?// °_o-Niveau 
iegen ; seine Multiplettkonstante betragt nach (30) etwa — 104 cm“. 
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* Die Formel fiir die Konstante C, in der Hut chissonschen Arbeit S. 415 oben ist nicht 
richtig ; sie muss lauten (im Falle von auf J normierten Sch wingungseigenfunktionen) : 


I i sat 
c= | eae 2(1 +a?) 
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O COCTOAHMH 477 MOJIEKYJIbI OF 


A. BYO u H. KOBAY 


Pe3rome 


MyJBTHIeTHoe pacmereHHe CNeKTpasbHOK JIMHHH MOJICKyJIbl Of, HaOmoAaemoe 
B cocTosHHn 4/7,, KoTOpoe, KaK @yHKHHA BpallaTeIbHOrO KBaHTOBOTO YHCIa OTMMYAeTCA 
oT POpMyNHI, BLIBeAeHHOM Ha OcHOBe Teopuu BaH-Oreka, HacTOmmeH paodoToH OObACHAeTCA 
BOSMyIeHHem cocTomHua 4/7 pBymMa cocTosHHsmn */Z,4TO WaeT XOpoulee CoBMaqeHHe C 9KCTIE- 
PHMCHTaIDHEIMH WaHHbiMH. H3 TeOopHH WaHHOrO BOSMyLeHHA MOKHO 3afaTb MpHOMMKeHHOE 
TOIOKKEHHE M MYJIbTHTIETHY!O KOHCTaHTy 40 cHx Mop HeusBecTHoro */I-Tepma. 


hyn 
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M. G. Mayer [1] explained the empirical nuclear spin values by assuming 


that of'two orbits in the same nucleus that of higher angular momentum has 


the higher (less negative) pairing energy. In her calculations to check this hypo- 
thesis, the potential between two nucleons being assumed in the shape of a 6 
function she obtained for the absolute value of the pairing energy belonging 
to the orbit determined by 


Prij a (27 i 1) In (1) 


where 


In = | Rai (r) r? dr (2) 
0 


and R,, (r) and r denote the normalized radial eigenfunction and the radius 
vector, respectively. 

(In what follows we shall neglect the factors independent of n, 1, j, as 
we are only interested in the relative value of the pairing energies belonging 
to the individual energy levels.) 

Thus, to determine the pairing energy we need to know the energy eigen- 
functions belonging to the individual quantum states. 

Hitherto the calculations have been carried out in two different ways. 
M.G. Mayer [1] estimated the value of the integral (2), whereas A. J. M. Hitel.- 
cock [3] used oscillator eigenfunctions; i.e. in his calculation the integral 
in the pairing energy actually depends only on the value of n, 1, but not on the 
value of j. In this way the increase with] of the pairing energy in the case of 
identical / is trivial. 

Thus a strong spin-orbit interaction seems a suitable assumption for 
our calculations, as in this case the energy eigenvalues and, at the same time, 
the eigenfunctions depend not only on n, | but also on j. Thus we obtain for the 


integral (2) : 
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oo 


Tay = [Ray (r) r? dr . (3) 


In our calculations we used G. Szamosi’s results [4]. Calculating the 
potential of the nucleus, which was assumed to be of constant density, from 
a Yukawa interaction potential the expression for the potential becomes 


ae 3gA [1—G tare SAT if r<R 
Vin) = iy 
3gA shuR |e“. 
oe eee Ge et, R L r>R 
Saal GeO re ea 


where A is the number of nucleons, and g the coupling constant of the meson 
field. 4 is related to the mass of the 2 meson in the well-known manner. In the 
Hartree approximation the one-particle Schrédinger equation belonging to the 
quantum number (n, I, 7) is of the form 


2M 
Ani + 95 9 


eee 
eee ev"| an. 


Assuming a spin-orbit interaction the potential becomes 


( 7 
| wie l iff jol4+— 
r dr 
V(r) = 
Ve (reece (+1) if es ee 
ip Ne 2 
where 
2 
“= ae 
Mec 


In the above expressions j and E,,,; denote the one-particle energy eigenfunction 
and the energy eigenvalue, respectively. M is the mass of the particle and 
A a numerical factor characteristic of the strength of the spin-orbit interaction. 


Separating angular and radial part in » the differential equation obtained for 
the function 


®,1; (r) = r Ray (r) 


can be solved numerically if g and A are given. 


ere 
on 
4 
4 
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Taking into account that the radial eigenfunctions in the integral are 
normalized we obtain for the absolute value of the pairing energy from (1) 
and (3) 


[O4,)(0) Far 
Bay 1) 


[| Gy @ arp 
0 


‘The calculations were carried out for the nucleus with Z = 60 and N = 83 
(Z number of protons, N number of neutrons), the parameter values being 
g = 5e and 4 = 10. The results of the calculation are listed in the Table. 


' . TABLE 
(The »competing« levels are marked by{). 


Level | A ei 
Tere 1,105 0,230 
2 Pa), 0,059 0,237 
ips: 0,061 0,122 

fled, 0,042 0,251 
\2 s, le 0,121 0,243 
Wed 0,045 0,180 
EE. 0,033 0,264 
2 Das. 0,058 0,233 
Ge Ate 0,038 0,229 
‘ Pi), 0,062 0,124. 
l gy), 0,029 0,289 
2 ds, 0,040 0,238 
Sy. | 0,126 0,251 
l‘g;,, | 0,033 0,264 
dye 0,043 0,172 
Phy 0,026 0,308 
oe. 0,030 0,243 
ths: 0,029 0,297 


The results show that the assumption for the pairing energy was justified : 
of orbits of the same azimuthal quantum number that of higher angular momen- 
tum has indeed the higher pairing energy. 

The same result, i.e. that it is really the orbit of higher angular mo- 
mentum which has a higher pairing energy, was obtained in the case of levels 
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which are, according to experiment «competing» [5]. Thus it is possible that 
the orbit of higher angular momentum lies energetically lower, namely 


Pray, > Pas, Prgy, > Papy,, Pre, > Peasy, 


There are only two exceptions to this rule, the levels 2p,, and lf... 
Calculating, however, the level order without taking into account the pairing 
energy, it was already found that these levels were interchanged as compared 
to the experimental results. Besides, the difference between the two pairing 
energies is so slight that it does not exceed the error limit of the calculation. 

These calculations are of mere informational character as in the derivation 
of (1) the interaction potential of the two nucleons was supposed to be a 6- 
function. It can be assumed, however, that the application of the Yukawa 
potential, apart from some constant factors, would yield essentially the same 
results. 

I take this opportunity to thank G. Szamosi for having called my atten- 
tion to this problem. 
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DER ENERGIE-IMPULS-TENSOR 
DER STRAHLUNG IN DIELEKTRIKA ~ 


Von 
G. MARX und K. NAGY 


INSTITUT FUR PHYSIK DER ROLAND EOTVOS UNIVERSITAT, BUDAPEST 


(Vorgelegt von K. F. Novobaétzky — Eingegangen: 29. VI. 1954) 


Die elektromagnetische Strahlung steht selbst in vollig durchsichtigen 


- Medien ‘in einer standigen Wechselwirkung mit der polarisierbaren Materie : 


das elektromagnetische Feld ruft eine wechselnde elektrische und magnetische 
Polarisation hervor, wahrend die Strahlung der entstehenden Dipole die Kigen- 
schaften der primar einfallenden Strahlung verindert. Hieraus folgt, dass die 
aus dem Vakuum einfallende Strahlungsenergie auch in Medien, die nicht 
standig Energie absorbieren, in einem gegebenen Augenblick nur teilweise 
in Form von elektromagnetischer Energie anwesend ist, wahrend sie teilweise 
als die potentielle und kinetische Energie der polarisierten (deformierten) 
Molekiile, also vom makroskopischen Gesichtspunkt als mechanische Energie 
erscheint. Mehrere Forscher wiesen darauf hin, dass der die Erhaltungssitze 
befriedigende, vollstandige Strahlungs-Energie-Impuls-Tensor S;, eben mit 
Ricksicht auf die obigen Ausfiihrungen im Inneren der polarisierbaren Materie 
vom Energie-Impuls-Tensor Tj, des elektromagnetischen Feldes zu unter- 
scheiden ist : zu Tj, muss man namlich jenen Teil t;, des mechanischen Energie- 
Impuls-Tensors addieren, der die vom elektromagnetischen Feld erregten 
Spannungen beschreibt. Das Ziel der vorliegenden Arbeit ist die Herleitung 
der durch die Feldstarken ausgedriickten kovarianten Form dieses Tensors. 

Mit der Bestimmung des vollstaindigen Strahlungstensors S;, haben sich 
bereits mehrere Forscher befasst. H. Ott [1] und F. Beck [2] bildeten den 
Strahlungstensor S;, auf Grund des Minkowskischen Energie-[mpuls-Tensors 
(um dessen Asymmetrie zu beseitigen), wahrend G. Marx und G. Gyérgyi [3 | 
vom Abrahamschen Energie-Impuls-Tensor ausgingen und auf Grund der 
ponderomotorischen Kriafte zu der auf ruhende Medien und ebene Wellen 
besch ankten Form von Sj; gelangten. K. Nagy [4] verallgemeinerte dann diese 
Form mit einer Lorentz-Transformation auch auf den Fall bewegter Dielektrika. 
Die Untersuchungen von Beck und Marx fiihrten trotz des verschiedenen Aus- 
gangspunktes zu einem im Wesen gleichen Ergebnis, was fiir die Richtigkeit 
ihrer Uberlegungen spricht. 

Bei der Herleitung des vollstandigen Strahlungs-Energie-Impuls-Tensors 
sei von dem Abrahamschen, Energie-Impuls-Tensor ausgegangen, der auch die 
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auf die im Dielektrikum entstehenden Polarisationsstréme wirkenden Lorentz- 
schen Krafte beriicksichtigt. Seine Form ist 


n?—1 


be 


T= Z| Fr Cr — EPs Gre — (Fi, Fr — Fr Fru) ue. (1) 
An 4 

Es werden, die iiblichen Bezeichnungen verwendet: F;, steht fir den aus den 
Vektoren © und % gebildeten Tensor, G;, fir den aus den Vektoren D und 9 
gebildeten Tensor, u; ist die in Lichtgeschwindigkeitseinheiten gemessene 
konstante Vierergeschwindigkeit des Dielektrikums, n der Brechungsindex des 
Mediums, “ die magnetische Permeabilitat des Mediums, s; die Viererstrom- 
dichte und F; = F;,u,. Die Maxwellschen Feldgleichungen lauten : 


9 Gix 


== Ah IES; 2 
= (2) 


SMe eer 0 Fix 


= 0. 3 
0 x; O xx Ox, ( ) 


Die die Feldgréssen miteinander verkniipfenden materiellen Gleichungen sind : 
1 
Ge Sa Ute eM rid) eae (4) 


Die nachstehenden Uberlegungen sollen sich zuniichst auf ruhende Dielektrika 
beschranken. Bezeichnet man die aus den raumartigen Komponenten des 
Tensors —Tj, gebildeten Maxwellschen Spannungen mit $, so gelangt man 
in einem homogenen, ladungs- und stromfreien (also durchsichtigen) Medium 
bekanntlich zu 


) 
Div F = — 
Ot 


n2 


_ Exg |. (5) 


4 c 


Man erhalt die Kraftdichte, wenn man auch die aus den Komponenten von 
Tia gebildete Feldimpulsdichte 


in Betracht zieht : 


: : oi n? — 1 
[= DiwT — = 
Ot a 47 ¢ Exo]. 
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Diese Kraft wird vom Feld auf das polarisierte Dielektrikum ausgeiibt. Wenn 
das Dielektrikum befestigt ist, dann vermag die ponderomotorische Kraft 
weder den Impuls des Dielektrikums zu erhéhen noch eine Arbeit an ihm zu 
verrichten. In diesem Falle muss also die im Dielektrikum erregte mechanische 
Spannung t die elektromagnetische Kraftwirkung gemiss dem Zusammenhange 


; t 


) [=e xs| 6) 


doce. 


kompensieren. Der mechanische Spannungstensor t hat also einerseits die 

Gleichung (6) zu befriedigen und andererseits natiirlich symmetrisch zu sein. 

Wenn man dies beriicksichtigt, dann kann man ganz einfach durch einen Ver- 

gleich der G1. (5) und (6) den Spannungstensor durch die Feldstarken ausdriicken: 
' 


Z. (7) 


Die Komponenten von t bilden hierbei die raumartigen Komponenten von —t;x. 
Die tibrigen Komponenten haben im Sinne der obigen Ausfiihrungen in einem 
starr befestigten Dielektrikum zu verschwinden. Zieht man die Verknipfung 
von & und Tj, sowie von t und t;, in Betracht, so lasst sich schreiben : 


n? — ] 
2 


tin = — Tins wenn oe ieeen ag phar te tg 0. 


n 
Diese fiir ein ruhendes Dielektrikum giiltige Relation kann zusammenfassend 
auch folgendermassen geschrieben werden : 


2 
tig = — (64, + wert) (Bas + tu te) Tre - (8) 
nN 


Der jetzt erhaltene Zusammenhang ist also fiir ruhende Medien giiltig, driickt 
aber daneben auch einen Tensorzusammenhang aus, so dass er in jedem Inertial- 
system Giiltigkeit besitzen muss. Durch die Gl. (8) wurde also die allgemein- 
giiltige kovariante Form von t;, erhalten. Hieraus lasst sich durch Heranziehung 


von 
Six = Tin + tix (9) 


auch der vollstandige Strahlungs-Energie-Impuls-Tensor gewinnen, u. zw. sowohl 
fir den Fall einer beliebigen elektromagnetischen Strahlung als auch fiir den 
Fall eines sich mit beliebiger Geschwindigkeit fortbewegenden Dielektrikums. 
Bei Anwendung der Gl. (1), (8) und (9) kann Six, auch als Funktion der 
Feldstarken ausgedriickt werden. Durch’ einfache Substitution ergibt sich 


- « a 
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: Tuy pea 

Sa =a [Fr Oe 2G) + ae [1 2a) a Bae Fr 


ib fi a 

EEE (n?—1) F, Fit Fr Fs \ : 
2 2 J 

(Durch Ubergang auf dreidimensionale Bezeichnungen kann bestatigt werden, 

dass dieser Zusammenhang mit der in [4] benutzten Form ibereinstimmt.) 

Die Spur des Tensors betragt nun — im Gegensatz zu der von T;, — nicht mehr 


Null, wodurch zum Ausdruck gebracht wird, dass ein Teil der Strahlungsenergie 
bereits als Energie einer tiber eine Ruhemasse verfiigenden Materie anwesend 


(10) 


ist: 


Lo SSS Srr 


a ae | 
— . : E Wopee oes Be] >0 : 
An un* | 4 2 


Bei Anwendung der Gl. (10) und der Feldgleichungen (2) und (3) kann auch 
unmittelbar durch Differenzieren bewiesen werden, dass S;, in einem durch- 
sichtigen, homogenen Medium divergenzfrei ist. An den inhomogenen und leit- 
fahigen Stellen des Dielektrikums (n const, ~~ const, s; ~ 0) weisen die 
Strahlungsenergie und der Impuls Senken auf : 


0 Six il 1 los 
aie En ee 1) up oe ee Pere = 
0 xK a Page ates 2a 82 n? : Xi 

1 1 
Fis hy Re 
82 n? uu? | 2 0 x; 


In diesem Falle haben aber die Abtrennung von t;, vom mechanischen Energie- 
Impuls-Tensor und infolgedessen auch der Strahlungstensor S;, bereits ihre 
Bedeutung verloren. Wo indessen die Energie nicht permanent absorbiert 
wird, sondern nur zwischen dem Feld und der polarisierten Materie oszilliert, 
kann S;, zur Beschreibung der dynamischen Eigenschaften der durch das 
Dielektrikum hindurchgehenden Strahlung (so z. B. der sich im Dielektrikum 
bewegenden Photone [4]) benutzt werden. 
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